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Preface 


This  thesis  was  stimulated  by  research  in  Tactical  Target 
Identification  performed  by  the  United  States  Air  Force  Systems 
Command's  Rome  Air  Development  Center  (RADC/OCTM),  Griff iss  AFB,  New 
York.  Engineers  there  are  investigating  the  possibility  of  iden¬ 
tifying  aircraft  from  two  dimensional  cross  range  versus  slant  range 
image  plots  of  the  discrete  electromagnetic  scatterers  on  the  air¬ 
craft.  To  identify  the  aircraft,  it  i..  thought  of  as  a  collection  of 
simple  geometrical  objects  whose  image  plots  are  known  at  different 
aspect  angles.  For  a  known  aircraft,  the  image  plots  of  the  indivi¬ 
dual  geometrical  objects  are  "pieced"  together  at  their  individual 
aspect  angles  with  their  spatial  position  observed,  with  respect  to 
some  reference  point  on  the  aircraft.  Thus,  what  is  produced  is  a 
composite  image  plot  of  the  discrete  point  scatterers  on  the  aircraft. 
The  assumption  made  above  is  that  secondary  multiple  scattering  is 
negligible  compared  to  primary  scattering. 

The  goals  of  this  thesis  are  to  obtain  a  practical  method  of 
obtaining  the  backscattered  signal  for  any  known  incident  signal  for 
the  perfectly-conducting  sphere  and  perfectly-conducting,  infinitely- 
thln  circular  disk  and  also  to  investigate  the  sensitivity  of  the 
backscattered  signal  to  changing  aspect  angle. 

I  ackowledge  with  thanks  the  advice  and  encouragement  received 
from  my  thesis  advisor,  August  Golden,  Jr.,  Capt.,  USAF  and  to  my 
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thesis  readers,  Joseph  W.  Carl,  Maj.,  USAF,  and  Professor  Raymond 
Potter  at  the  Air  Force  Institute  of  Technology.  Also,  my  thanks  and 
appreciation  go  to  my  sponsor,  Mr.  Richard  Wood  of  the  Surveillance 
Division,  Rome  Air  Development  Center  and  to  Drs.  D.  L.  Moffat  and 
D.  B.  Hodge  of  the  Ohio  State  University.  Dr.  Hodge's  computer 
simulation  for  the  scattering  from  the  circular  disk  was  extremely 
helpful.  Also,  Capt.  Golden's  many  helpful  suggestions  were  inva¬ 
luable  to  the  completion  of  this  research. 

Finally,  I  would  like  to  thank  my  wife,  Kathy,  whose  patience 
and  understanding  were  of  immeasurable  value  throughout  my  AFIT 
program. 
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Abstract 

The  problem  of  identifying  the  perfectly  conducting  sphere 
and  perfectly-conducting,  inf initely-thin  circular  disk  via  the 
backscattered  far  field  from  the  object  with  a  known  incident  signal 
is  examined.  A  method  utilizing  Synthetic  Aperture  -adar  principles 
to  obtain  slant  range  versus  cross  range  image  plots  of  the  discrete 
point  scatterers  on  an  object  is  discussed.  A  computer  program  to 
calculate  the  far  field  backscattering  from  the  sphere  is  developed. 
Another  computer  program,  found  in  the  literature,  which  calculates 
the  far  field  backscattering  from  the  inf initely-thin,  circular  disk, 
is  extended  for  ka>15.  In  the  above  programs,  the  backscattering  is 
obtained  as  a  function  of  frequency  and  treated  as  the  frequency 
response  of  a  linear  system.  The  input  to  this  system  is  five  cycles 
of  10GHz  R-F.  The  return  signal  is  the  output  of  the  system  and  is 
obtained  by  using  conventional  linear  system  theory.  Since  a  sphere 
is  symmetric,  the  same  return  is  obtained  at  all  aspect  angles.  This 
return  from  the  sphere  is  examined.  However,  returns  from  the  cir¬ 
cular  disk  are  different  at  different  aspect  angles.  These  returns 
from  the  disk  are  examined. 
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PREDICTED  MICROWAVE  ELECTRO-MAGNETIC  BACKSCATTERING 
RETURNS  FOR  SIMPLE  REFLECTIVE  TARGETS 

I.  Introduction 

Background 

With  the  use  of  Synthetic  Aperture  Radar  processing,  it  is 
theoretically  possible  to  obtain  two-dimensional  slant  range  versus 
cross  range  image  plots  of  the  locations  of  the  discrete  point  scat¬ 
tered  on  simple  geometrical  shapes.  These  plots,  it  is  hoped,  will 
show  a  rough  outline  of  the  object  and  will  aid  in  its  identification. 
As  was  stated  in  the  Preface,  an  aircraft  will  be  thought  of  as  a 
collection  of  these  simple  geometrical  objects  whose  image  plots  are 
known  at  different  aspect  angles. 

The  two  simple  geometrical  shapes  considered  here  are  the 
sphere  and  the  circular  disk.  The  sphere  is  unique  in  the  fact  that 
at  any  aspect  angle  the  return  signal  will  always  be  the  same.  This 
is  due  to  spherical  symmetry.  Thus,  at  any  specific  aspect  angle,  the 
image  plot  obtained  from  the  sphere  will  be  the  same  as  the  image  plot 
obtained  at  any  other  aspect  angle.  As  will  be  shown  in  Chapter  II, 
to  obtain  cross  range  information,  the  return  signal  must  change  with 
changing  aspect  angle.  Thus,  for  the  sphere,  only  slant  range  infor¬ 
mation  can  be  obtained.  The  circular  disk,  however,  has  the  advantage 
of  returning  an  aspect-dependent  signal.  Therefore,  cross  range 
Information  can  be  obtained  for  the  disk.  The  circular  disk  also 
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possesses  a  known  backscattered  frequency  response  which  is  treated  as 
the  frequency  response  of  a  linear  system.  The  backscattered  signal 
is  the  output  of  the  system  and  the  incident  signal  is  the  input  to 
the  system.  Although  the  backscattered  signal  from  the  sphere  provi¬ 
des  no  cross  range  information,  its  backscattered  signal  can  be 
obtained  in  the  same  way  as  the  backscattered  signal  from  the  disk. 

The  method  used  for  identification  of  the  unknown  shape  is  to 
compare  the  image  obtained  from  the  unknown  shape  to  previously 
recorded  images  of  different  known  shapes  at  known  aspect  angles. 
Hopefully,  one  of  the  known  image  plots  will  match  up  with  the  unknown 
plot.  Taken  that  this  is  the  case,  the  unknown  object  will  be  iden¬ 
tified  as  the  object  which  corresponds  to  the  known  matching  image 
plot. 


Problem  and  Scope 


The  specific  problem  studied  in  this  thesis  is  how  many 
distinct  backscatter  signature  returns  are  obtained  over  a  specified 
angle  variation  around  the  circular  disk.  The  sphere  is  spherically 
symmetric  and  thus  will  backscatter  the  same  return  signal  for  all 
incident  aspect  angles.  For  this  reason,  the  section  on  the  sphere 
was  used  to  show  how  the  return  signal  was  obtained  by  knowledge  of 
the  frequency  response  of  the  scattering.  The  return  signal  from  the 
disk  is  obtained  in  the  same  manner.  In  this  case,  however,  the  fre¬ 
quency  response  of  the  scattering  Is  aspect  dependent. 

Only  the  sphere  and  the  circular  disk  are  considered  in  depth 
since  they  are  the  only  simple  geometrical  shapes  for  which  the 
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frequency  response  of  the  far  field  scattering  can  be  numerically 
calculated. 

General  Approach  and 
Presentation 

The  general  approach  used  in  the  research  was  to  firstly 
understand  how  two-dimensional  images  of  the  discrete  point  scat¬ 
tered  on  simple  geometrical  objects  can  be  obtained.  Synthetic  Aper¬ 
ture  Radar  processing  was  seen  to  aid  in  the  development  of  a  slant 
range  versus  cross  range  image  of  the  discrete  point  scatterers  of  the 
object. 

Secondly,  there  had  to  be  discovered  if  there  were  any 
geometrical  shapes  for  which  the  backscattered  field  as  a  function  of 
frequency  was  numerically  calculable.  Once  the  sphere  and  the  cir¬ 
cular  disk  were  found  to  possess  this  quality,  they  were  the  only  sha¬ 
pes  concentrated  on.  As  was  alluded  to  in  the  background  section, 
this  backscattered  field,  as  a  function  of  frequency,  was  treated  as 
the  frequency  response  of  a  linear  system.  The  incident  signal  is  the 
input  to  the  system  and  the  backscattered  signal  is  the  output.  For 
the  sphere,  only  the  backscattered  signal  at  one  aspect  angle  is  exa¬ 
mined  since  the  sphere  is  symmetric.  For  the  disk,  however,  the 
returns  from  different  aspect  angles  were  examined. 

The  point  of  obtaining  the  backscattered  signals  from  the 
sphere  and  the  disk  is  that,  theoretically,  they  can  now  become  the 
inputs  to  the  processor  to  produce  the  two-dimensional  plots. 
Unfortunately,  because  of  time  constraints,  this  could  not  be 
attempted. 
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The  sequence  of  presentation  of  the  main  chapters  in  the 

sis  is: 


I.  Introduction 

II.  Signal  Processing-Two  Dimensional  Images 

III.  Scattering  from  the  Perfectly  Conducting 
Sphere  in  the  Far  Field 

IV.  Scattering  from  the  Perfectly  Conducting, 

Infinitely-Thin  Circular  Disk  in  the  Far  Field 

V.  Conclusions,  Findings  and  Recommendations. 
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II.  Signal  Processing  Two-Dimensional  Images 

The  processing  used  to  obtain  two-dimensional  images  of  the 
scattering  centers  of  objects  in  the  far  field  is  essentially 
Synthetic  Aperture  Radar  processing  (Ref  3,  Ref  8).  The  following 
is  an  explanation  of  the  processing  used  to  obtain  a  two-dimensional 
image  of  two  discrete  point  scatterers  in  the  far  field  (Ref  7). 

Figure  1  is  representative  of  the  situation  of  two  discrete 
point  scatterers  in  the  far  field.  Scatterer  1  is  at  point  one  and 
scatterer  2  is  at  point  two.  Imagine,  however,  that  instead  of  the 
scatterers  moving  around  the  radar,  that  the  radar  moves  around  the 
scatterers.  The  radar  is  at  point  three.  With  no  loss  in  generality, 
let  scatterer  1  be  located  at  zero  cross  range  and  let  scatterer  2  be 
located  at  cross  range  R^.  Also  let  the  distance  from  the  radar 
to  scatterer  1  be  constant  and  be  denoted  Rg.  Let  the  distance  from 
the  radar  to  scatterer  2  be  varying  and  be  denoted  R.  Also  let 
A9/2  and  0  denote  the  angles  shown  in  Figure  1.  The  angle  A9/2s0° 
and  9s90°.  Therefore,  the  cross  range  axis  is  perpendicular  to  the 
line  from  point  one  to  point  three. 

The  two  dimensional  image  that  will  be  obtained  will  be  a 
plot  of  slant  range  versus  cross  range.  Slant  range  can  be  obtained 
easily  by  noting  the  time  delay  of  the  signal  from  the  radar  to  each 
point  scatterer  and  back  to  the  radar.  Therefore,  if  the  time 
delay  to  the  i^  scatterer  and  back  is  denoted  t^,  the  slant  range 

Rsi  ’  1s 


5 


(1) 


where 

c-  speed  of  light  in  free  space 

As  can  be  seen,  the  processing  required  to  obtain  slant  range 
is  quite  simple  in  principle.  On  the  other  hand,  the  processing 
required  to  obtain  cross  range  is  more  complicated.  The  cross  range 
of  scatterer  2  will  be  examined  since  the  processing  involved  to 
obtain  its  cross  range  is  representative  of  the  processing  involved 
to  obtain  the  cross  range  of  any  discrete  point  scatterer  on  the 
target. 

Assuming  the  transmitted  signal  from  the  radar  is 
sT(t)  .  eJ“c‘  (2) 

the  signal  received  back  from  scatterer  2  will  be 

sR2(t)  ■  ej(a)ct'2kRs2)  (3) 

where 

k  =  f  <4> 

Also,  assuming  that  cross  range  values  will  be  much  smaller  than 
slant  range  values,  can  be  approximated  by 

r$2  s  ro"Rc2  cos  9  (5) 
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therefore , 


sR2(t)  *  c>cl  •  e-J2kR0  •  eJ2Rc2  cos  0 


(6) 


note  that 


fl-a  -  JU+Z  +  M 
0  “  9(t)  Rq  2  2 


(7) 


therefore 


cos  e  =  cos  (f  +  •y-  -  -^-t) 


-  -  sin  [~2  -  -j^t) 

However,  if  -  -^-t)  is  small,  then 
L  K0 


(8) 


q  -  v*  A0 
COS  6  s  -  -g- 


(9) 


Therefore, 


W‘>  ■ e 


jo)  t  -j2kRg  -jkA©Rc2  J  R 

S  A  a  A  a  A  a  0 


2kR~v 
j-^-t 


(10) 


At  the  receiver,  there  will  be  a  bank  of  matched  filters. 
Each  of  them  will  be  matched  to  a  specific  cross  range  on  either 
side  of  zero  cross  range.  There  will  also  be  a  filter  matched  to 
zero  cross  range. 

The  bank  of  matched  filters  is  represented  in  Figure  2 
where  Ax  is  the  cross  range  resolution  and  T  is  the  dwell  time. 

If  there  Is  a  discrete  scatterer  at  a  cross  range  which  is  an 
integer  multiple  of  Ax,  the  outputs  of  the  filters  matched  to  those 
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integer  multiples  will  roll  off  with  a  sin  x/x  distribution.  For 


example,  suppose  there  is  a  scatterer  at  cross  range  R  .  Therefore 

v 


2kR  v 

jo)  t  -j2kRn  -jkA6R  j— g^-t 

sR2<t)  =  e  c  •  e  0  .  «  c  •  e  R0 


(ID 


The  e^c*  term  can  be  suppressed  without  loss  in  generality.  There¬ 
fore, 


sR2(t)  =  e 


-j2kRr 


-jkAQR, 


2kR  v 
j  b~  t 


(12) 


The  output  of  the  filter  matched  to  Rc2  would  be 


2kR/.v  2kRr9v 

T/2  -jk(2R  +A6R  )  dt 

/  e  •  e  0  •  e  0 

-T/2 


.  e-^<2vieRc>  dt 

-T/2 


(13) 


=  e 


-jk(2RQ+A0Rc) 


T  sinCf^(Rc"Rc2)T:i 

t^(Rc-Rc2>TJ 


A  plot  of  the  magnitude  of  the  above  equation  as  a  function  of  R 

V 

is  shown  in  Figure  3.  If  Rc=  R^,  a  maximum  of  T  will  be  the  out¬ 
put  of  the  filter  matched  to  rc2.  If  RC?‘RC2»  the  output  of  that 
same  filter  will  fall  off  as  a  sinc(x)  distribution  as  a  function 


To  determine  cross  range  resolution  the  criteria  that  will 

be  used  is  that  if  the  output  of  the  filter  matched  to  Rc2  is  above 

0.707T  in  magnitude,  then  there  is  a  scatterer  at  cross  range  R^: 

sinc(x)=0.707  at  approximately  1.4  Rad.  Therefore,  AR  =R  -  R  0  must 

c  c  c  c 

be  found  in  Eq.  (14)  and  multiplied  by  two  to  obtain  Ax. 


sinCf^RJ] 


=  0.707 


Therefore 


^  «eT  s  iA 


AR.  = 


1.4XRfl  0.22XRQ 


c  2-rrvT  vT 

0.44XRn 


Ax  =  2AR.  = 


c  VT 


(14) 


(15) 


(16) 


(17) 


After  slant  range  and  cross  range  of  each  scatterer  are 
obtained,  a  two  dimensional  plot  is  made  for  a  single  target  at 
a  single  aspect  angle  by  locating  these  scatterers  on  a  display. 
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III.  Scattering  from  the  Perfectly  Conducting 
Sphere  in  the  Far  Field 

The  electromagnetic  backscattering  from  an  object  can  be 
thought  of  as  a  linear  system  characterized  by  a  frequency  response 
which  is  the  normalized  backscattered  E-field  as  a  function  of  fre¬ 
quency  (Ref  2:6).  Therefore,  in  theory,  the  frequency  spectrum 
of  the  backscattered  signal  can  be  obtained  by  multiplying  the 
Fourier  Transform  of  the  incident  signal  by  the  frequency  response 
of  the  system.  The  backscattered  signal  is  the  inverse  Fourier 
Transform  of  its  frequency  spectrum.  Thus,  by  conventional  linear 
system  theory,  the  scattered  signal  from  an  object  can  be  obtained. 

The  reason  for  obtaining  the  backscattered  signal  from  the 
object  is  that  it  can  now  be  inputed  to  the  processor  described  in 
Chapter  II  in  order  to  obtain  a  two-dimensional,  slant  range  vs.  cross 
range  image  of  the  discrete  point  scatterers  on  the  object.  The 
process  is  represented  in  a  block  diagram  in  Figure  4. 

Frequency  Response  of  the  Sphere. 

The  Mie  series  solution  of  the  scattering  from  the  sphere  at 
a  point  p  as  a  function  of  ta  is  (Ref  14:142): 

Es(p,w)  =  EQ  ?  (AnMoln+BnNeln)  (18) 

n=l 

where. 
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Roln  =  sTn“e  hJ(kr)p^C0S  0)(cos  <f>)e 


-hj(kr)^pj(c°s  9)(sin  4»)4> 


Nein  =  ^(kr)pj(cos  0)(cos 


+  FWtkrhS(kr)]lP!(cos  9){cos  *)9 


FritHl  dlRT  Ckrhl(kr)jpl{c°s  0)(cos  +)* 


pj(x)  -  Associated  Legendre  Polynomial  of  degree  n  and  order  1. 


”Jf —  ~(x2-l)n 

2n*n!  dxn 


h*(x)  -  Spherical  Hankel  Function  of  degree  n  and  order  1. 


=  0_(x)  +  jY  (x) 

n  n 
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(23) 


1 


Y_(x)  =  JnM  cos  ^  *  J_nM 
n  sin  (mr) 


n  not  an  integer 


A 


n 


-(-J) 


n  2n+l 
n(n+l ) 


3n(ka) 

h’(ka) 


(24) 


,  .  d  [kajn(ka)l 

n  ^  t  i \H+1  2n+l  d(ka) 

Bn  nTn+T)  d  r.  .1, 

d^kahn<ka>J 


(25) 


However  .assuming  r»a,  e“ju)t  time  dependence  and  using  far 
field  approximations,  with  "a"  the  radius  of  the  sphere  (Ref  14:143): 


Ene 


£s(p,«)  =  -° 


jkr 


where 


P —  [cos  <f>S-j (e)0  -  sin  0S2(0)<f>] 


S1<0)  -  ?  +  jB  d  p'( 

n=l  sin  e  n  de  n 


cos  e)] 


S2(8)  •  (-j)n+,CAn  ^'(cos 


Pn(cos  0) 


0 


] 


let, 


(26) 


(27) 


(28) 


F(0,<J>)  T  =  cos  4>S1  (e)e  -  sin  <f>S2(@)e 


(29) 


In  the  backscattering  direction,  approximations  are  made 
for  F(0)  for  four  regions:  the  low-frequency  region,  the  lower 
resonance  region,  the  upper  resonance  region  and  the  high 
frequency  region  (Ref  14:146-150): 
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In  the  low  frequency  region  (ka<0.4)  the  approximation  for 
F(0)  is  (Ref  14:150): 

F(0)=f(ka)2  (30) 

therefore 

E  eJ’kr 

rEs(a>)  =  •  f(ka)3  (31) 

In  the  lower  resonance  region  (0.4<ka<0.8)  the  approxima¬ 
tion  for  F(0)  is 


F(0)  =  |(ka)3[l-^(ka)2^(ka)4 


6,651,923 

11,907,000 


(ka)6] 


+  4(ka)6[l  +  fCka)2] 


‘  j 
e 


(32) 


where 


x  =  |(ka)3[l-^(ka)2+^(ka)4 


6,651,923 

11,907,000 


(ka)6] 


(33) 


(34) 


y  -  {(ka)6[l  +  f(ka)2] 
therefore 


rEs(oj) 


E(jix2+y2  j(kr  +  tan-1  *£) 
- -  .e 


(35) 


In  the  upper  resonance  region  (0.8<ka<20)  the  approximation 
for  F(0)  consists  of  a  specular  component,  F°(0),  added  to  a  creeping 
wave  term,  Fc(0).  Therefore, 

F(0)  *  F°(0)  +  Fc(0)  (36) 

The  approximations  of  F  (0)  and  F  (0)  are: 
o  ka  "J2ka  i 

F°(0)  =  -f  e  Cl-2kV  1  (37) 

_  4  J^(ka  -  !■) 

Fc(0)  =  (kar  e  6  [1.357588  + 

_  2 

+  (0.741196  +  jl.283788)(ka)3]  • 
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1 

•  exp  [-(ka)3(2. 200002  -jl. 270172) 

1 

+  (ka)  3(0. 445396  +  jO. 257150)] 

-2 

+[0.695864  +  (0.964654  +  j 1 . 670829 )( ka )  3]  • 

1 

•  exp  [-<ka)3(7. 014224  -  j4. 049663) 

-  (ka)  3(0. 444477  +  jO. 256619)] 

2 

-  [0.807104  +  (0.798821  +  jl.383598)(ka)  3]  • 

1 

•  exp  [-(ka)^( 5. 048956  -  j2. 915016) 

-  (ka)  3(0. 312321  +  jO. 180319)]  (38) 


let 

F°(0)  +  Fc(0)  =  Vx2+y2  etan_1  x  (39) 

where 


x  =  Re{F°(0)  +  Fc(0) } 


(40) 
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y  =  Im(F°(0)  +  Fc{0)} 


(41) 


therefore 


rEs(u>)  =  Eq 


j(kr+tan-1  J) 
e  * 


(42) 


In  the  high  frequency  region  (ka>20)  the  approximation  for 

F(0)  is 


F(0)  =  -  ^kae‘j2ka  (43) 

therefore 

rEs(w)  =  -  e-J[kr-2kal 

(44) 

.  .  lol  gj[kr-2ka] 


A  computer  program  is  listed  in  Appendix  A  which  calculates 
the  frequency  response  of  the  far  field  E-field  backscattering  from 
a  perfectly  conducting  sphere.  To  give  an  idea  of  the  shape  of 
the  backscattered  E-field  frequency  response,  a  sample  of  the 
response  is  shown  in  Figure  5  as  a  function  of  ka  up  to  ka=30. 
Since  the  scattering  from  the  sphere  is  aspect  independent,  the 
response  shown  in  Figure  5  is  valid  at  any  aspect  angle. 


-  H 

I  ~  II  'ar»'*r»mrtl 
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The  facility  for  obtaining  Fourier  Transforms  and  Inverse 
Fourier  Transforms  in  this  thesis  is  the  Fast  Fourier  Transform  (FFT). 
An  explanation  as  well  as  a  listing  of  the  computer  program  for  the 
FFT  is  provided  in  Appendix  C.  A  review  of  the  proper  selection  of 
FFT  parameters  is  provided  in  Appendix  D. 

Process  to  Obtain  Backscattered  Signal 

The  incident  signal  used  is  a  pulse  of  five  cycles  of  10GHz 
R-F.  The  selection  of  the  proper  way  to  sample  the  incident  pulse  is 
explained  in  Appendix  D.  From  Appendix  D,  the  pulse  is  shown  to  be 
sampled  ten  times  each  cycle,  or  fifty  times  each  pulse.  The  pulse  is 
shown  in  Figure  6  centered  at  t=0.  The  pulse  is  shown  split  in  two 
parts.  One  part  is  at  t=0  while  the  other  part  is  at  t=l0.24  nsec. 

This  is  done  because  the  FFT  cannot  see  negative  time.  The  FFT  assumes 
periodicity  of  the  time  signal  and  its  frequency  spectrum.  Therefore, 
the  pulse  would  repeat  around  t=10.24  nsec  (NOTE:  1024  samples, 
samples  are  0.01  nsec  apart).  An  expanded  view  of  the  left  half 
of  the  pulse  is  shown  in  Figure  7.  The  spectrum  of  the  pulse  is  shown 
in  Figure  8.  Note  that  the  spectrum  is  simply  a  sinc(x)  function 
shifted  in  frequency  and  centered  on  10GHz.  As  outlined  in  Appendix  D, 
the  frequency  response  of  any  object  must  go  out  to  the  250th  sample, 
with  Af  between  samples  being  97, 656,000Hz,  to  fully  enclose  the 
appreciable  frequency  components  of  the  incident  10GHz  pulse.  The 
frequency  response  for  the  sphere  is  shown  in  Figure  9,  with  the 
left  half  frequency  reversed,  conjugated  and  shifted  up  in  fre¬ 
quency.  The  conjugation  occurs,  because  when  this  frequency  response 
is  inverse  transformed,  the  result  should  be  a  real  impulse  response. 
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Figure  6.  Incident  Pulse  of  R-F  at  10  GHz 
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It  will  be  real  if  the  frequency  response  has  even  magnitude  and  odd 
phase.  The  real  part  of  the  inverse  Fourier  Transform  of  the  fre¬ 
quency  response  in  Figure  9,  is  shown  in  Figure  10. 

In  Figure  11  is  shown  the  real  part  of  the  result  that  is 
obtained  when  the  frequency  response  of  the  sphere  is  multiplied  by 
the  Fourier  Transform  of  the  10GHz  incident  pulse  and  inverse  Fourier 
Transformed.  Figures  12-14  show  the  same  process  done  on  a  video  pulse 
of  the  same  width  as  the  10GHz  pulse.  Figure  12  is  the  video  pulse 
with  Figure  13  its  spectrum,  and  Figure  14  the  real  part  of  the  back- 
scattered  signal.  Figures  11  and  14  will  be  compared.  Because 
the  video  pulse  has  frequency  components  down  very  low  in  frequency 
where  the  frequency  response  of  the  sphere  is  oscillatory,  the 
response  obtained  by  linearly  convolving  this  pulse  with  the  impulse 
response  of  the  sphere  is  not  exactly  the  video  pulse  returned  as 
Is  shown  in  Figure  14.  However,  the  appreciable  frequency  components 
of the  frequency  spectrum  of  the  10GHz  pulse  are  out  in  frequency 
where  the  frequency  response  of  the  sphere  is  essentially  flat. 
Therefore,  as  shown  in  Figure  11,  the  10GHz  pulse  is  returned  almost 
intact  with  only  a  time  delay.  The  slight  negative  response  at  the 
end  of  the  graph  in  Figure  11  also  holds  much  significance.  This 
is  the  return  from  the  incident  signal  creeping  around  the  sphere. 

The  significant  point  is  that  as  the  incident  pulse  gets  shorter 
and  shorter  in  duration,  more  and  more  physical  features  of  the 
object  that  the  pulse  is  illuminating  will  become  noticeable.  In 
the  case  of  the  sphere,  there  are  only  two  returns;  the  specular 
return  and  the  creeping  wave  return.  An  unfortunate  result  is 
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Figure  10.  Real  Part  of  Inverse  Fourier  Transform  of  the 
Frequency  Response  of  the  Scattering  from  the 
Sphere  28 
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Figure  13.  Spectrum  of  Video  Pulse 


31 


0.00  5.00  10.00  15.00  20.00  25.00  30. 

_ FREQUENCY  F _  *109  (Hertz) 


REFLECTED 


0.00  2.00  4.00  6.00  8.00  10.00  12.00  14.00 

_ TIME'  T _ *10~9  (seconds) _ 


that  as  the  frequency  of  the  pulse  is  increased,  the  creeping 
wave  term  will  become  less  apparent.  So  only  the  specular  return 
will  be  noticeable.  At  lower  frequencies,  the  creeping  wave 
return  will  become  more  apparent.  Refer  back  to  Iq.  (17)  in 
Chapter  II,  written  here  for  convenience. 

0.44XRn 

Ax  =  2ARC  = - ^  (17) 

As  the  frequency  is  decreased  in  order  to  notice  the  creeping  wave 
return,  the  wave  length,  X,  will  be  increased  and  thus  increase 
Ax  and  thus  decrease  the  cross  range  resolution.  As  the  frequency 
is  increased  in  order  to  decrease  X  and  Ax  and  therefore  increase 
the  cross  range  resolution  of  the  discrete  point  scatterers,  the 
creeping  wave  return  will  disappear.  Thus,  for  all  practical  purposes, 
only  the  specular  return  will  be  detected. 


IV.  Scattering  from  the  Perfectly  Conducting, 

Infinite! y-Thin  Circular  Disk  in  the 
Far  Field 

As  was  stated  in  the  Introduction,  the  sphere,  until  recently 
was  the  only  geometrical  object  for  which  numerical  scattering  results 
could  be  obtained.  Although  the  rigorous  eigenfunction  solution  of 
the  scattering  from  the  thin  circular  disk  has  been  available  in  the 
literature  for  some  time  (Ref  1),  only  limited  numerical  results 
have  been  available  due  to  the  complexity  of  the  numerical  computations 
involved.  However,  in  February  1979,  a  report  out  of  the  Electrosciences 
Laboratory  at  the  Ohio  State  University  was  published  for  the 
Department  of  the  Navy,  Office  of  Naval  Research,  which  made  avail¬ 
able  the  rigorous  eigenfunction  solution  of  the  plane  wave  scatter¬ 
ing  problem  from  the  infinitely-thin  circular  disk  (Ref  9).  The 
author.  Dr.  D.  B.  Hodge,  develops  a  computer  program  to  obtain  the 
far  field  scattering  from  the  thin  metallic  disk.  The  program  allows 
an  incident  plane  wave  of  arbitrary  incidence  angle  and  polariza¬ 
tion  and  calculates  the  far  field  radar  cross  section  as  well  as 
the  scattered  E-field.  The  geometry  of  the  disk  is  shown  in  Fig¬ 
ure  15.  The  direction  of  the  incident  wave  is  specified  by  k1  and 
9q.  The  incident  direction  is  restricted  to  the  x-z  plane.  Because 
of  circular  symmetry,  this  is  no  restriction.  The  angle  "a"  is 
the  polarization  angle.  When  a=0°,  the  incident  field  is  parallel 
to  the  x-z  plane.  When  a=90°,  the  incident  field  is  perpendicular 
to  the  x-z  plane.  The  radius  of  the  disk  is  "a."  While  the  incident 


Figure  15.  Geometry  of  Disk 
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field  is  restricted  to  the  x-z  plane,  the  scattered  field  can  be 
in  any  direction  specified  by  8s  and  $  .  However,  since  backscat- 
tering  is  only  considered,  0s=0g  and  4>s=0.  A  listing  of  Dr.  Hodge's 
program  is  contained  in  Appendix  B  along  with  a  sample  I/O  and  an 
explanation  of  how  to  input  values  to  the  program. 

The  program  is  valid  up  to  ka=15.  A  high-frequency  approxi¬ 
mation  was  used  to  extend  the  program  past  ka=15.  This  approxima¬ 
tion  is  a  modification  of  the  Sommerfield-MacDonald  technique  by 
Ufimitsev  (Refs  16  and  17).  For  O°<0g<9O°  and  0s=6q,  the  theta  com¬ 
ponent  of  the  E-field  normalized  by  the  radius  of  the  disk  and 
divided  by  two  is  (Refl4:514-515). 


-3 


A^(9)  ( 


1-sin  6n  r,  1+sin  6n 
2  sin  +  A2  ^®^2  sin  0A 


•  JjUka  sin  0q)  - 
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(0)  ( 


1-sin  0n 

- 

2  sin  0Q; 


(45) 


9  1+sin  0n 


Jg^ka  sin  0g) 


where 


S  (n)  =  /0  Sin  -|t2  dt 

°||  =  47n"2!E0l2 


(46) 


(47) 


(48) 


(49) 


To  limit  the  extent  of  the  programs  listed  in  Appendix  B, 
ct=0,  in  Figure  15.  Since,  such  is  the  case,  only  Eg  need  be  con¬ 
sidered  since  will  be  zero.  The  program  used  to  implement  the 
above  equations  for  Eg  and  aj  j  is  in  Appendix  B  with  an  output 
following  the  first  program.  The  computer  program  that  incorpo¬ 
rates  the  extension  to  the  first  program  is  the  third  program 
listed  in  Appendix  B. 

Sensitivity  of  Returns  from  the  Circular 
disk  at  Five  Different  Aspect  Angles' 

The  following  discussion  is  divided  into  three  sections. 

The  first  section  concerns  backscattering  from  the  circular  disk  at 
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1°  off  normal  incidence.  The  second  section  is  a  sensitivity  study 


of  the  backscattering  from  the  disk  at  71°,  72°,  73°,  74°,  and  75° 
off  normal  incidence.  The  last  section  examines  backscattering 
from  the  disk  at  89°  off  normal  incidence. 

Backscattering  from  the  Circular  Disk  at  1°  Off  Normal. 

Figure  16  shows  the  frequency  response  of  the  backscattered  E-field 
from  the  circular  disk  at  1°  off  normal.  The  large  oscillatory  portion 
of  the  response  is  the  backscattering  in  accordance  with  the  physical 
optics  model.  If  the  frequency  of  the  incident  E-field  is  constant 
and  is  directed  at  a  specific  aspect  angle  off  the  normal  of  the 
disk,  a  lobing  pattern,  quite  similar  to  an  antenna  lobing  pattern, 
will  be  created  with  the  maximum  of  the  major  lobe  pointed  at  the 
same  angle  off  normal  as  the  incident  E-field  and  180°  around  the 
disk  from  the  incidence  direction.  This  pattern  is  the  magnitude  of 
the  scattered  E-field  as  a  function  of  aspect  angle.  As  the  fre¬ 
quency  is  increased  and  the  incidence  angle  is  kept  constant,  the 
beam  widths  of  the  lobes  of  the  lobing  pattern  become  smaller  and 
consequently  maxima  and  minima  of  the  lobing  pattern  pass  by  the 
fixed  incident  angle  resulting  in  the  oscillatory  nature  of  the  fre¬ 
quency  response  of  the  backscattered  E-field.  The  oscillatory 
response  riding  on  the  physical  optics  response  is  due  to  diffraction 
off  the  edges  of  the  disk.  If  this  response  were  to  be  plotted  out 
by  itself,  it  would  be  seen  to  be  oscillatory  in  nature  although  of 
much  smaller  magnitude  than  the  physical  optics  response.  Figure  17 
is  a  plot  of  the  phase  of  the  frequency  response  whereas  Figure  16 
is  a  plot  of  the  magnitude  of  the  frequency  response.  Figure  18 
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is  a  plot  of  the  backscattered  signal  at  1°  off  normal  incidence 
when  a  video  pulse  0.5  nsec  wide  with  a  pulse  repetition  frequency 
of  100MHz  is  incident.  This  "impulse  response"  is  seen  to  have  a 
doublet  at  the  leading  edge  of  the  response.  This  contribution  is 
due  to  a  term  proportional  to  ka  in  the  scattered  electric  field  and 
is  predicted  by  both  physical  optics  and  the  geometrical  theory  of 
diffraction  (Ref  10:559).  Although  the  response  in  Figure  18  is 
not  truly  the  impulse  response  since  the  incident  pulse  is  not  an 
impulse  but  a  video  pulse  of  width  0.5  nsec,  the  response  obtained 
closely  resembles  results  obtained  in  the  literature  (Ref  10 : p . 559 ) . 
The  backscattered  signals  obtained  when  a  pulse  of  five  cycles  of 
10GHz  R-F  is  incident  at  1°  is  shown  in  Figure  19.  The  effect  of  the 
doublet  is  still  seen  as  the  phase  of  the  signal  immediately  to  the 
right  of  t=0  has  a  phase  of  180°  relative  to  the  incident  signal 
shown  in  Figure  6. 

Backscatterinq  from  the  Circular  Disk  from  71°,  72°,  73°,  74°, 
and  75°  off  Normal .  As  was  stated  in  the  Introduction,  the  specific 
problem  studied  in  this  thesis  is  how  many  distinct  backscatter 
signature  returns  are  obtained  over  a  specified  angle  variation 
around  the  circular  disk.  The  vehicle  for  accomplishing  this  was  to 
find  the  angle  for  which  the  backscatter  frequency  response  had  the 
lowest  minimum  at  10GHz  (remember  the  incident  pulse  is  five  cycles 
of  10GHz  R-F).  The  angle  is  73°  and  was  obtained  with  the  use  of  the 
combination  program  or  3rd  program  in  Appendix  B.  The  angle  variation, 
mentioned  above,  was  centered  on  73°.  This  variation  in  angle  was 
from  71°  to  75°  in  one  degree  increments.  The  frequency  responses 
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which  correspond  to  the  five  angles  are  shown  in  Figures  20-24. 

The  backscattered  signals  at  these  angles  are  shown  in  Figures  25-29. 
If  more  time  were  available,^  backscattered  signals  would  have 
been  inputed  to  the  signal  processor  outlined  in  Chapter  II,  but 
because  of  time  constraints  only  the  backscattered  frequency  responses 
and  backscattered  signals  will  be  compared.  Their  similarities  and 
dissimilarities  will  be  pointed  out  as  well  as  conclusions  that  can 
be  dervied  from  the  figures. 

In  Figure  20  is  shown  the  frequency  response  of  the  back- 
scattering  when  the  incidence  direction  is  71°  off  normal.  Note  the 
minimum  is  at  approximately  8GHz.  As  theangle  is  increased  from 
71°  through  75°,  this  minimum  will  move  to  the  right  and  would  pass 
through  10GHz  when  the  incident  angle  equals  73°  to  approximately 
13GHz  when  the  angle  equals  75°. 

If  images  were  able  to  be  made  from  the  backscattered  signals 
corresponding  to  angle  variations  around  the  five  angles  discussed, 
one  of  two  things  would  be  true.  If  the  images  were  virtually  the 
same,  then  only  one  of  the  images  would  be  recorded.  Thus,  this 
would  reduce  the  number  of  image  plots  to  characterize  the  object. 

If  the  image  plots  were  quite  dissimilar,  then  all  the  image  plots 
obtained  would  definitely  be  needed  and  the  angle  variations  would 
further  have  to  be  subdivided  to  see  if  the  image  plots  corres¬ 
ponding  to  the  smaller  angle  variations  were  similar  or  not.  This 
process  is  repeated  until  image  plots  which  were  virtually  the 
same  were  obtained.  At  this  point,  only  the  previous  image  plots 


are  retained  for  comparison  to  the  image  plot  obtained  from  the 
unknown  object. 

In  Figures  25  to  29  are  shown  the  backscattered  signals 
when  the  angle  of  incidence  is  71°,  72°,  73°,  74°  and  75°  off 
normal,  respectively.  Notice  that  the  backscattered  signal  when 
0q=73°  is  the  smallest  return  out  of  the  5  returns.  This  makes  sense 
since  at  73°  the  frequency  response  of  the  backscattering  has  a 
minimum  at  10GHz  which  is  the  frequency  of  the  incident  signal. 

As  0q  is  increased  and  decreased  around  73°  the  backscattered  signals 
are  shown  to  increase  in  amplitude.  This  indicates  that  a  valley  in 
the  frequency  response  has  been  passed  through.  At  the  right 
side  of  the  curves  shown  in  Figures  25-29  there  is  another 
response.  This  is  conjectured  to  be  the  diffracted  return  with 
the  dominant  curve  on  the  left  being  the  specular  return. 


Backscattering  from  the  Circular  Disk  at  89°  off  normal. 

In  Figure  30  is  shown  the  frequency  response  when  00=89°.  This 
case  is  considered  to  illustrate  what  the  return  signal  looks  like 
at  almost  edge-on  incidence. 

In  Figures  31  and  32  are  shown  the  backscattered  signals 
at  0q=89o  for  the  10GHz  R-F  Pulse  and  video  pulse.  Note  that  the 
amplitudes  of  the  returns  are  reduced  from  the  amplitudes  of  the 
returns  centered  around  0q=73°  by  2  orders  of  magnitude.  It  can  thus 
be  concluded  that  the  diffracted  returns  are  always  smaller  than  the 
physical  optics  returns.  It  is  interesting  to  note  the  curve  on 
the  right  side  of  Figure  32.  This  return  quite  possibly  could  be 
the  pulse  traversing  the  disk  and  being  diffracted  back.  From 
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Figure  20.  Frequency  Response  of  the  Scattering 
from  the  Circular  Disk  with  0q®71° 
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Figure  21.  Frequency  Response  of  the  Scattering 
from  the  Circular  Disk  with  8q=72° 
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Figure  22.  Frequency  Response  of  the  Scattering 
from  the  Circular  Disk  with  8g=73° 
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Figure  25.  Reflected  Signal  from  the  Circular 
Disk  with  R-F  Pulse  Incident  at  71° 


Figure  26.  Reflected  Signal  from  the  Circular 
Disk  with  R-F  Pulse  Incident  at  72° 
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Figure  30.  Frequency  Response  of  the  Scattering 
from  the  Circular  Disk  with  0O=89° 
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Figure  31.  Reflected  Signal  from  the  Circular 
Disk  with  R-F  Pulse  Incident  at  89° 
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Figure  32.  Reflected  Signal  from  the  Circular 

Disk  with  Video  Pulse  Incident  at  89° 
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Figure  32  it  can  be  observed  that  the  time  delay  between  the  beginning 
of  the  first  return  and  the  beginning  of  the  second  return  is 
approximately  5.20  nsec.  It  can  be  calculated  that  the  distance 
in  space  that  corresponds  to  a  time  delay  of  5.2  nsec  is  equal  to 
approximately  1.56  m.  This  means  that  the  incident  pulse  has  to 
travel  an  additional  1.56  m  to  produce  the  second  return  than  it 
had  to  travel  to  produce  the  first  return.  The  first  return  in 
Figure  32  is  considered  the  physical  optics  return  from  the  edge 
of  the  disk.  The  second  return  in  Figure  32  is  conjectured  to  be 
the  return  caused  by  the  incident  pulse  traversing  the  disk  and  then 
traveling  half  way  around  the  disk  to  the  point  of  incidence. 

With  the  radius  of  the  disk  equal  to  0.30  m,  the  additional  distance 
that  the  second  return  travels  is  approximately  1.54  m  which  provides 
a  time  delay  of  approximately  5.20  nsec  which  corresponds  to  the 
time  delay  between  returns  in  Figure  32. 
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VI.  Conclusions.  Findings,  and  Recomnendations 


Conclusions 

There  are  two  major  conclusions  that  can  be  derived  from  this 

thesis. 

The  first  conclusion  is  that  no  cross-range  information  can  be 
obtained  from  the  perfectly-conducting  sphere  due  to  the  aspect  inde¬ 
pendent  nature  of  the  backscattered  return.  The  second  conclusion  is 
that,  with  the  10GHz  pluse  incident,  due  to  cross-range  resolution 
constraints,  only  the  specular  scatterer  and  not  both  the  specular  and 
creeping  wave  scatterers  will  be  detected  on  an  image  plot.  Thus,  the 
image  plots  obtained  from  the  processing  discussed  in  Chapter  II  would 
only  show  one  scatterer  and  would  not  aid  in  the  identification  of  the 
sphere. 

Findings 

The  backscattered  signals  centered  around  =  73°  for  the 
circular  disk  are  quite  similar.  The  above  indicates  that  maybe  the 
image  plots  that  would  be  obtained  from  angle  variations  around  71°, 
72°,  73°,  74°  and  75°  off  normal,  quite  possibly  could  be  quite  simi¬ 
lar.  If  this  were  the  case,  the  five  image  plots  could  be  replaced 
by  one  image  plot.  If  such  were  not  the  case,  then  these  five  image 
plots  would  definitely  be  needed  in  the  series  of  image  plots  used  to 
characterize  the  circular  disk. 
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Recommendations 

It  is  recommended  that  the  processor  described  in  Chapter  II 
be  implemented  and  that  the  number  of  distinct  image  plots  that  are 
obtained  for  the  circular  disk  be  investigated.  The  effect  of 
increasing  or  decreasing  the  frequency  or,  equally,  decreasing  or 
increasing  the  wavelength  of  the  incident  E-field,  on  the  cross  range 
resolution  in  reference  to  the  disk  should  be  investigated.  In  other 
words,  are  the  diffraction  returns  off  the  edges  of  the  circular  disk 
attenuated  as  the  frequency  is  increased  and  is  Ax,  as  defined  in 
Chapter  II,  too  large  to  resolve  the  diffraction  scatterers  as  the 
frequency  is  decreased. 
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Appendix  A 


Explanation  and  Program  Listing  for  the 
Frequency  Response  of  the  Sphere 

The  following  listing  is  the  computer  program  that  was  used  to 
generate  Figures'  6,  7,  8,  9,  10,  and  11.  The  program  is  a  straight¬ 
forward  application  of  the  equations  for  the  different  frequency 
regions  as  outlined  in  Chapter  III.  However,  a  few  words  will  be 
stated  to  explain  the  plotting. 

First  of  all,  the  plotting  was  done  using  a  library  which  is 
unique  to  the  computer  system  at  the  Air  Force  Institute  of 
Technology.  The  library  of  subroutines  used  is  called  SAPL.  The 
subroutine  within  SAPL  used  to  do  the  plotting  is  called  HGRAPH.  As 
an  example,  generation  of  one  of  the  graphs  will  be  explained. 

Figure  10  is  a  graph  of  the  real  part  of  the  impulse  response 
of  the  sphere.  The  computer  program  statements  which  generate  this 
graph  are  statements  1140  to  1220  inclusive.  The  array  IE  which  has 
17  locations  is  used  for  labeling.  IE(1)  to  IE(8)  are  used  to  label 
the  upper  righthand  corner  box.  IE(9)  and  IE(10)  are  used  to  label 
the  abscissa  and  IE(ll)  and  IE( 12)  are  used  to  label  the  ordinate. 

IE( 13)  to  IE( 17 )  are  used  to  title  the  graph  (NOTE:  10  hollerith 
characters  fit  in  one  array  location).  This  title  goes  on  the  outside 
of  the  graph  box  along  the  ordinate  coordinate.  The  CALL  PLOT  state¬ 
ments  merely  position  the  pen  on  the  plotter.  The  CALL  HGRAPH  state¬ 
ment  calls  the  correct  subroutine  for  plotting.  The  argument  list 
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r 


contains  seven  arguments.  In  the  example,  V  Is  the  array  of  x-axis 
values  of  dimension  1026.  ROESR  is  the  array  of  y-axis  values  of 
dimension  1026.  Both  arrays  are  two  more  locations  than  1024  for 
Internal  scaling  purposes.  The  "1024"  is  the  number  of  points  to  be 
plotted.  The  IE  is  the  labeling  array.  The  "1"  specifies  internal 
scaling  by  the  program.  The  first  "0“  means  connect  successive  points 
by  a  straight  line.  And,  finally,  the  second  "0"  means  to  draw  dots 
at  points  as  opposed  to  other  symbols  (i.e.,  triangles,  circles, 
etc.). 

The  preceding  is  all  that  is  needed  to  plot  graphs  using  the 


HGRAPH  subroutine  in  the  SAPL  library. 


I  0  0=  PROGRAM  I  MAGE  1  £  <■  I NPUT »  Dir  PUT .  PLOT) 

I I  0=  PERL  C 1  ,  C£ «  E 1 ,  E£ » F 1  »  F£,  G 1 «  G£ ,  H 1 » H£ , I1,I£,J1>  J£ » K 1 » Ve‘ » 

1  £  0=  8L  1  ,  L£  ,  M 1  *  M£  ,  N 1 » N£ »  P 1  *  P£  t  Kfi 

130=  INTEGER  Z 

1 4 0=  COMPLEX  C»E»F»G»H»I»J»K*L*M»N»P 

1 5  0=  8 ,  PULSE < 1 0£4> * PULSER < 1 0£4> » ROES  < 1 0£4> 

160=  D I MENS I ON  X < 1 0£6>  » Y < 1 0£6> » U < 1 0£6> ,V<1 0£6> *  ID <  1 7> 

170=  Sf  I A  <17>  » IB  <17>  » IC<17>  »  IE  <17>  » ROESR  (1  0£6>  ,  POES  I  <10£6> 
130=  S» IF  <17> 

1 9  0=C 

£00=C  FIX  DIRMETER  OF  SPHERE » ft ,  AND  DISTANCE  TO  TARGETED. 

£1  0=C 

220=  •  R=0.  3 

£30=  R0=0. 0 

240=  B=0. 45 

£50=C 

£6 0=C  COMPUTE  MHRT  THE  SCATTERED  FIELD!.  ROES-,  fii  R  FUNCTION 

£70=C  OF  WAVE  NUMBER'.  IS  FOP  THE  LOU  FREQUENCY  REGION,  THE 

£8 0=C  LOWER  RESDNPNCE  REGION,  THE  UPPER  RESONANCE  REGION, 

£90=C  AND  THE  HIGH  FREQUENCY  REGION. 

3  0  0=C 

310=  FREQ=97656  000 

3£0=  DO  10  Z=1 , £51 

33 0=  KR=  <  < 0 .  3>  ♦  < £>  ♦  O.  141 59>  ♦FREQ'*- <Z- 1 )  >  ✓  <3 0 0 0 0 0 0 0 0> 

340=  IF  <KA. EQ. 0> GO  TO  55 

350=  C1  =  0.0 

36  0=  C£= <KR^RD> ✓ <A> 

37  n=  C=CMPLX < C 1 » C£> 

380=  E 1  =  <3 .  s£ .  >  ♦  <KA^3 .  >  ♦  <  1 .  -  <5 .  /54 .  )  ♦  <KR^£ .  > 

390=  8  +  <17,  /900.  >  ♦  <k'A^4.  >  -  <6651 9£3.  ^11 907000.  >♦  <KA^*6.  V 

400=  E£=  <  1  ./'£•>  ♦  <KR^6  .  >  ♦  <  1 .  +  <6.  ^5.  >  ♦<KR^£.  >  > 

410=  E=CMPLX  <E 1 9 E£> 

4£0=  F 1 = 1 . 0 

430=  F£=l .  /  <£. ♦KA> 

44  0=  F=CMPLX <F 1 , F£> 

450=  G1  =  0.  0 

460=  G£=- <£. ♦KA> 

47  0=  G=CMPLX  <6 1 ,  G£0 

430=  H1  =  0. 0 

490=  H£=  <3.  141593>  ♦  <KA-  <1 .  •'‘6.  >  -1 .  ) 

5  o  0=  H=CMPLX  <H 1 » H£> 

510=  11  =  1.357588  +  <0.  741 196>  ♦  <KR~  <-£.  s2.  >  > 

5£0=  I£= <1 . £83788> ♦ <KR*» <-£. s2. > > 

530=  I=CMPLX  <1 1 , I£> 

54  0=  J 1  =-  <KA*^  <  1 .  s'3  .>>♦<£.£  0  0  0  0£>  + 

55  0=  $  <KR^  <- 1 .  s’ 3 .  >  >  ♦  <  0 . 445396> 

560=  J£= <Kft»^ <1 . s3. > > ♦ <1 . £7017£>  + 

570=  8  <KA~<-1 .  /3:  >  >  ♦  <0.  £57150> 

580=  J=CMPLX <J1 , J£> 

590=  K  1  =  0. 695864  ♦  <0.  964654)  ♦  <KR^  <-£..  s3.  >  > 

600=  K£=—  <1 . 6708£9>  ♦  <KR^<-£.  s’3.  >  > 

610=  K=CMPLX  <K 1 , k£> 

6£0=  Ll=- <KA^  <1 . s3, >  > ♦(T. 014££4>  - 

630=  8  <Kfl^*  <-l .  -'3.  >  >  ♦  <0.  444477) 

64 0=  L£=  <KA#^  < 1 . S3 .  >  >  ♦  <4 . 049663>  - 

650=  8  <KA*^  <-l .  ^'3.  >  >  ♦  <0.  £566 19> 


66  0= 
670= 
680= 
68  0= 
700= 
710= 
720= 
730= 
740= 
750= 
760= 
770= 
780= 
790= 
800= 
810= 
820= 
830= 
84  0= 
850= 
86  0= 
870= 
880= 
890= 
90  0= 
910= 
920= 
930= 
94  0= 
9^  0= 
960= 
970= 
980= 
990= 
1000= 
1010= 
1020= 
1030= 
104  0= 
1050= 
1060= 
1070= 
1080= 
1090= 
1100= 
1110= 
1120= 
1130= 
1140= 
1150= 
1160= 
1170= 
1180= 
1190= 
1200= 


L=CMPLX CL 1 « L2> 

M 1  =  0 . 8  07 104  +  C  0 . 79882 1  >  ♦  CKR^  C -2 .  /3 .  >  > 

M2=  c  1 . 383598}  ♦  CKfi-*^  c-£.  /3.  >  > 

M=CMPLX CM 1 , M2> 

N 1  =-  C KR-**  <  1 .  /3 .  >  >  ♦  C5 .  048956}  - 
8  CKfi^  <- 1 .  .'-'3 .  }  }  ♦  <  0.312321} 

N2= CKR^ < 1 . ^3 .>}♦ <2 .915016}  - 
$  CKfi^  c- 1 .  -"3 .  }  }  ♦  <  0.  1 8 03 1 9} 

N=CMPLXCN1»N2> 

P1  =  0.  0 

P2=  C  CKfi^RO)  ✓  <fi}  }  -  <2.  }  ♦  CKfi} 

P=CMPLX  CPI • P2} 

IF  CKfi.  LE.  0.4}  GO  TD  20 
IF  C KR .  6T .  2 U .  U.‘*  GO  TD  5 0 
IF  CKfi.  LE.  0.8>60  TD  30 
I F  CKfi . GE . 0 . 9} GD  TD  4  0 
2 0  RDES  CZ}  =  CB-'-Kfi}  ♦  CKfi-^3 . }  ♦  CCEXP  CC}  } 

GD  TD  10 

30  ROES  CZ}  =  C  C  C2.  /3.  >.  ♦£)  /  CKfi}  >  ♦  CCEXP  CC}  >  ♦  CE} 

GD  TD  10 

4 0  ROES  CZ}  =  C  C  C  C2 .  /3 . }  ♦B)  CKfi}  }  ♦  CCEXP  CC}  }>♦ 

8  C-  CKfl''2 .  }  ♦  CCEXP  CG}  }  ♦F  +  CKfi**  C4 .  /3  .}>♦ 

$  CCEXP  CH}  >♦  Cl }♦  CCEXP  CJ}}  +  <K)  ♦  CCEXP  CL>  >  - 
3  CM}  ♦ CCEXP CM} } } 

GD  TD  10 

50  ROES  CZ}  =  C-B-z'3. }  ♦  CCEXP  CP}  } 

GD  TD  10 

55  RDES  CZ}  =  <0. »  0. } 

10  CONTINUE 

HD  60  Z=1 * £51 
R=PERL CRDES CZ} > 

B=RIMRG  CF'DES  CZ}  } 

RDES CZ} =CMPLX Cfi » -B} 

X  CZ} =RERL CRDES  CZ} } 

Y  CZ} =R I MRG CRDES  CZ> } 

6  0  W  CZ}  =  C  C  C  CX  CZ} } ♦♦£ . }  +  C CY  CZ} } ♦♦£ .}}}♦♦ C 1 .✓£. } 

HD  65  Z=£5£?  774 
RDES  CZ}  =  CO. > 0 • } 

65  W  CZ}  =  0.  0 

BO  68  Z=1 >  £50 

RDES C 1 025-Z} =CMPLX CX  CZ+ 1 } r - CY CZ+ 1}}} 

68  CONTINUE 

CRLL  FFT  CRDES » 1 0» -1 . } 

DD  1  Z=1 f 1 024 

RDESP  cz:>  =F'ERL  CRDES  CZ}  > 

RDES I CZ} =fi I MRG  CRDES  CZ}  } 

VCZ}  =  CZ-1.}/'C10^U.} 

1  CONTINUE 


BRTfl  IE  Cl }■••'£ OH  PERL  PRPT  DF  / 

BfiTP  IE  C3} "20H  IMPULSE  RESPONSE 
BRTR  IE  C5} z£0H  DF  SPHERE  / 

IiflTfi  I E  C7}  ■''£  OH  BIBRRTDLDMED  ✓ 

DfiTR  IE  C9}  z£0H  T  ✓ 


BftTft  IE  Cl 1} ^2 OH  IMPULSE  RESPONSE  ✓ 
HfiTfi  IE  03}  /5  OH 
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1210=  CALL  PLOT  CO. * -4. . -3. >  $  CALL  PLOT  CO. * . 03* -3. > 

1  £3  0=  CALL  HGPAPH  1 V  •  ROE  -  F-  *  1  034  *  I E  *  1  *  0  *  0  * 

1330=  DATA  I F  '1  >  s 3  0 H  IMAGINARY  F'AF'T  OF  / 

1340=  DATA  IF  <'3>  ••'3 OH  IMPULSE  RESPONSE  ✓ 

1350=  DATA  I FC5  OH  OF  SPHERE  / 

1360=  DATA  IF  C?>,-£ OH  D I BARTOLOMEO  ✓ 

1370=  DATA  IFC9V-30H  T  ✓ 

1380=  DATA  I F  C 1 1 > OH  IMPULSE  RESPONSE  ✓ 

1 28  0=  DATA  I F  C 1 3>  s% OH 

/ 

1300=  CALL  PLOT  C 0 .  *  -4 .  >  -3 .  >  $  CALL  PLOT  CO.  *.  03*  -3.  > 

131 0=  CALL  HGPAPH CV* ROESI *  1 034* IF* 1  *  0* 0) 

1330=  CALL  FFT CEDES *  1 0* 1 . ) 

1330=  DO  69  Z=775» 1 034 

1 34  0=  X  CZ>  =REAL  CROES'  CZ)  ) 

1 35  0=  V <Z> =A I MAG CROES <Z> > 

1 36  0=  W  CZ)  —  C  C  C  CX  cz:>  )  ♦♦2 .  >  +  C  <Y  CZ)  )  ♦♦2 .  )  )  )  ♦♦  Cl.  /£ .  ) 

1370=  69  CONTINUE 

1380=  DO  70  Z=1 * 20 

1 39  0=  PR  I  NT ♦ ,  X  CZ>  >  V  CZ>  >  W  CZ> 

1400=  70  CONTINUE 

1410=  DO  80  Z=1 *  1 034 

1420=  V CZ> =FREG-*Z 

1430=  80  CONTINUE 

144  0=  DATA  I D  C 1  )  /£ OH  FREQUENCY  RESPONSE  ✓ 

1450=  DATA  ID  C3>  -"2 OH  OF  SPHERE  / 

1460=  DATA  ID  C5>  -X20H  D I  BARTOLOMEO  / 

1470=  DATA  I D  C7> /£ OH  ✓ 

1480=  DATA  I D  C9) /£ OH  F  ✓ 

1490=  DATA  I D  C 1 1 ) /£ OH  ROES CF)  ✓ 

1500=  DATA  ID  Cl 3) ^  5 OH  FREQUENCY  RESPONSE  DF  SPHERE 

✓ 

1510=  CALL  PLOT  CO.  *  -4.  *  -3.  >  $  -CALL  PLOT  CO.  *  .  03*  -3.  > 

1 53  0=  CALL  HGPAPH CV *  W *  1 034  *  I D  > 1 »  0 * 0> 

1530=  DO  90  Z=1 *  25 

1 54  0=  X= CZ+S5) ✓ < 1 <)♦♦ 11.) 

1 55 0=  PULSE  1  =S IN  C£^3.  14159*  Cl  1  0 .  >  ♦X) 

1560=  PULSES- 0.0 

1 57  0=  PULSE CZ) -CMPLX CPULSE 1  *  PULSES) 

1580=  90  CONTINUE 

1590=  DO  100  Z=36* 999 

1600=  PULSE 1=0.0 

1610=  PUL  SEE- 0.0 

1 63  0=  PULSE  CZ> -CMPLX CPULSE 1  *  PULSE£> 

1630=  100  CONTINUE 

1 64  0=  DO  1 05  Z= 1 0 0 0 *  1 034 

1 65  0=  X* CZ-999 - ✓ C 1 0*+ 1 1 . > 

1 66  0=  PULSE  1 =S I N  <3*3 . 1 4 1 59»  c 1 1 0 . > ♦X) 

1670=  PULSE2=0. 0 

1680=  PULSE <Z) -CMPLX CPULSE 1* PULSES) 

1690=  105  CONTINUE 

1700=  DO  110  Z=1 *  1 034 

1710=  X CZ)  =PEAL CPULSE  CZ>  ) 

1 72  0=  V  CZ)  =  CZ- 1 .  >  ✓  <  1  0^  1 1 .  ) 

1730=  110  CONTINUE 

1740=  DATA  I A  C 1  >  /•  2  OH  INCIDENT  / 

1750-  DATA  I A  C3> ^3 OH  PULSE  ✓ 
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176  0= 
1770= 
1780= 
1780= 
1800= 
/ 

1810= 

1880= 

1 83  0= 

184  0= 
1850= 
186  0= 
1870= 
1880= 
1890= 
1900= 
1910= 
1920= 
1930= 
194  0= 
1950= 
1960= 
1970= 

/ 

1980= 

1990= 

£000= 

£010= 

£0£0= 

£030= 

£040= 

£050= 

£060= 

£070= 

£080= 

£090= 

£100= 

£110= 

£1£0= 

£130= 

£140= 

£150= 

2160= 

/ 

2170= 

2180= 

£190= 

£200= 

£210= 

2220= 

£230= 

£240= 

2250= 

2260= 

2270= 

2280= 

2290= 

2300= 


BhTh  19  <  5  > -  2 OH  Dlf  9RT0L0UE0  / 

DR  T9  19  *:'7>  •••'20H  / 

D9T9  1 9  •:  9>  -'  2  OH  T  / 

D9T9  19  -l  1>  .-80H  PULSE  CT>  ✓ 

D9T9  19  a 3>  . '  5 OH  INCIDENT  PULSE  OF  PF  9T  10  GH 

C9LL  PLOT  C 0 . ? -4 . ? -3 . >  8  CRLL  PLOT  C  0 . ? . 03 * -3 . > 

C9LL  HGR9PH  <V?  X ?  1  024 ?  19?  1 «  0?  0.' 

C9LL  PLOT  <  0 . ? -4 . ?  -3 .  >  $  CRLL  PLOT  < 0 .  »  .  03  ?  -3 .  > 

CRLL  HGP9PH  *  V ?  X?  8  00?  19?  1  ?  0?  CO 
CRLL  FFT (PULSE? 10? l.> 

DO  120  Z=1 ? 1 024 
X  <Z>  =PE9L  < PULSE  CZ>  > 

V  <Z> =91 M9G (PULSE CZ>  > 

W  <2>  =  <  C  <  <X  <Z>  >  ♦♦£ .  >  +  <  CY  <z>  >  ♦♦£  .>>>♦♦ C 1 ./£. > 


120  V CZ>  =FPEQ^Z 

D9T9  IE(1>.-20H  TP9NSF0PM  OF  ✓ 

D9T9  IE  •' 3 >  -'20H  PULSE  ✓ 

D9T9  I E  C5>  ••"£  OH  D I BARTOLOMEO  / 

D9T9  I E  <7> /£  OH  ✓ 

D9T9  IE  C9>  •'20H  F  / 

D9T9  I B  C 1 1  >  •''£  OH  PULSE <F>  ✓ 


D9T9  I E  C 1 3> /5 OH  FOURIER  TR9NSFDPM  OF  PULSE 


140 


150 


C9LL  PLOT  C  0 . ? -4 . ? -3 . >  $  CRLL  PLDT  <  0 . ? . 03  ? -3 . > 
CRLL  HGP9PH C V ?  W ?  £ 0 0  ? IE? 1 ? 0 ? 0) 


DO  140  Z= 1 ? 1 024 

PULSER <Z>  =  CROES CZ> > ♦ CPULSE  <Z> > 

CONTINUE 

CRLL  FFT  CPULSER  ?  1  0  ?  - 1 .  > 

DO  150  Z= 1 ? 1 024 
X <Z> =RERL (PULSER  <Z> > 

V  CZ  ) =9 1  NRG ( PULSER  CZ>  > 


W  <z>  =  C  C  <  C X  <Z>  >♦♦£.>  +  <  < Y  CZ>  >  ♦♦£ . 
v  cz>  =  <:z- 1 .  >  ✓  c  i  i  .  > 

CONTINUE 

DRTR  ICa>/£0H  REFLECTED 
DRTR  I C  <:3>  2  OH  SIGNRL 
DRTR  IC  (5>  ''20H  DIERRTOLONEO 


>  >  ♦♦Cl 

/ 

/ 

/ 


-'£.  > 


DRTR  I C  C7>  ••'  £  OH 
DRTR  ICC9.:'  .'£0H 
DRTR  IC  Cl 1>  •■■£ OH 
DRTR  ICC13>/50H 


/ 

T  x 

PULSER  <T> 

REFLECTED  PULSE  OF  RF 


CALL  PLDT  CO.  ? -4.  ? -3.  >  S  CRLL  PLOT  CO. ?. 03? -3. > 
CRLL  HGRRF’H  (  V  ?  X  ?  1  0£4?  I C  ?  1  ?  0  ?  0> 

CRLL  PLOT  CO.? -4.  > -3 .  >  $  CALL  PLDT  CO.?. 03  ? -3 . > 
CRLL  HGR9PH  CV  ?  V  ?  1  024  ?  I C ?  1  ?  0 ?  CO 
CALL  FFT (PULSE ? 1 0 ? - 1 . > 

DO  160  Z=1 ? 1 024 
X  CZ> =RE9L (PULSE  CZ>  > 

1 6  o  v  <z>  =  cz- 1  :> "  c  i  o*t  1 1 .  :> 

CRLL  PLOT  CO.  ? -4 .  ? -3 .  >  8  CALL  PLOT C 0. ? . 03? -3. > 
CALL  HGRRPH CV?  X? 1 024? I R? 1 ?  0? 0> 

CALL  PLOTECN!) 

STOP 

END 

SUBROUTINE  FFTCX?M?XI> 
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2310=  COMPLEX  X < 1 > * U » W * T 

2320=  H=2»»M 

2330=  MV2=N'£ 


234  0= 

NMl=h— 1 

2350= 

J=1 

236  0= 

DO  7  1  =  1)  MM  1 

2370= 

I F  < 1 . 6E .  J>  GO  TO  5 

2^80= 

T=X  < J> 

2390= 

X  <  J>  =X  <  I  > 

2400= 

X  <I>  =T 

2410= 

5 

K=NV£ 

2420= 

6 

IF(K.6E.J>  GO  TO  7 

2430= 

J=.J-K 

244  0= 

K=K'2 

2450= 

GO  TD  6 

£460= 

7 

J=J+K 

£470= 

PI =3. 14159265358979 

£480= 

DO  £0  L= 1 ?  M 

£490= 

LE=£~L 

2500= 

LEl=LE/2 

2510= 

U=  <1.0)0.  0> 

2520= 

W=CEXP  (CMP LX  <  0 .  )  -X I  ♦P I  ■'-LE 1  >  > 

253  0= 

DO  20  J= 1 ) LE 1 

£540= 

DO  10  I=J) N) LE 

2550= 

IP=I+LE1 

2560= 

t=x  ( i  p:> 

2570= 

x  <  i  p:>  =x  <  i  >  -t 

2580= 

10 

X  <I>  =X  <I>  +T 

£590= 

£0 

1  l=IJ*l.l 

2600= 

IFCXI.GT. 0.>  RETURN 

2610= 

DO  30  1=1 *N 

£620= 

30 

X  <  I>  =X  <  I >  'N 

2630= 

RETURN 

Appendix  B 


Explanation  and  Program  Listings  for  the  Frequency 
Response  of  the  Far  Field  Scattering  of  the 
Perfectly-Conducting  Infinitely-Thin, 

Circular  Disk 

The  first  computer  program  shown  on  the  following  pages  cal¬ 
culates  the  far  field  radar  cross  section  and  associated  E-field  for 
the  infinitely  thin  circular  disk  (Ref  9).  A  sample  of  the  program 
I/O  is  shown  after  the  program  listing.  The  user  need  only  input 
eight  variables  for  the  initial  case  to  be  executed.  They  are: 

1.  KA  =  electrical  radius  of  the  disk  =  (2-rr/A) a . 

2.  THETA  INCIDENT  =  eQ  (degrees) 

3.  POLARIZATION  =  a  (degrees) 

4.  THETA  SCATTERED  =  0$  (degrees) 

5.  PHI  SCATTERED  =  a  (degrees) 

WHICH  VARIABLE  IS  TO  BE  INCREMENTED?  : 

if  1  through  5:  the  variable  associated  with  that  index 

above  will  be  incremented 

if  1  through  5:  then  only  the  initial  case  will  be 

calculated.  In  this  event,  the  remaining 
parameters  are  not  requested. 

TYPE  NUMBER  OF  CASES:  enter  the  total  number  of  computa¬ 
tions  to  be  performed. 
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WHAT  IS  THE  INCREMENT?:  enter  the  increment  by  which  the 
variable  selected  should  be  increased  after  each 
execution. 

All  inputs  may  be  made  in  free  format.  The  program  labels 
the  first  column  with  the  name  of  the  variable  to  be  incremented. 

The  second  and  third  columns  contain  the  cross  sections  associated 
with  the  e  and  <$>  components  of  the  scattered  field.  The  last  four 
columns  list  the  magnitudes  and  phases  (in  degrees)  of  both  the  e 
and  <p  components  of  the  normalized  scattered  electric  field. 

On  completion  of  the  cases  desired,  the  program  will  ask  for 
a  new  KA  and  proceed  as  before.  If  a  0  is  entered,  the  program 
will  terminate.  If  a  -1  is  entered,  a  brief  statement  of  the  problem 
geometry  will  be  given.  Following  this  statement,  the  program  will 
ask  for  a  new  value  of  KA. 

The  program  works  fine  up  to  KA=15.  For  KA>15,  a  high 
frequency  approximation  is  used  as  is  outlined  in  Chapter  IV.  The 
program  to  implement  the  high  frequency  approximation  is  called 
SCATTER  and  follows  the  I/O  of  the  first  program.  The  output 
for  0q=45°  follows  SCATTER'S  program  listing.  It  is  important  to 
note  that  in  relation  to  the  input  parameters  of  the  first  program, 
SCATTER  is  restrictive  in  the  sense  that  it  returns  only  the  theta 
components  of  cross  section  and  E-field  and  assumes  the  phi  compo¬ 
nents  of  cross  section  and  E-field  are  zero  (a=0).  Also  the  incident 
angle  is  equal  to  the  scattered  angle  (backscattering)  and  is 
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restricted  between  0°  and  90°.  It  is  also  restrictive  because  it 
will  only  increment  KA. 

Although  the  first  program's  output  and  the  output  from 
SCATTER  have  values  for  the  cross  section  and  E-field  from  KA=l-25, 
the  first  program  is  invalid  for  KA>15  and  SCATTER  is  invalid  for  the 
lower  values  of  KA.  What  is  needed  is  a  program  that  combines  the 
two.  This  combination  program  follows  the  output  from  SCATTER.  The 
output  from  the  combination  program  follows  its  listing.  Note  that 
for  KA<15,  it  uses  values  from  the  first  program  and  for  KA>15,  it 
uses  values  from  SCATTER. 

In  the  SCATTER  program,  the  Fresnel  Integrals  were  calcu¬ 
lated  using  the  trapezoidal  rule  for  approximation.  The  values 
obtained  from  the  FRESNEL  subroutine  agree  with  the  literature  out  to 
three  decimal  places  (Ref  11:34). 
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100* 
110= 
120= 
130= 
14  0= 
150= 
160= 
170= 
180= 
190= 
200= 
£10= 
£20= 
£30= 
£4  0= 
£50= 
£6  0= 
£70= 
£80= 
£90= 
300= 
310= 

32  0= 

33  0= 

34  0= 
350= 
360= 
370= 
380= 
390= 
400= 
4 1  0= 
420= 
430= 
44  0= 
450= 
460= 
470= 
480= 
490= 
500= 
510= 
520= 
530= 
54  0= 
550= 
560= 
570= 
580= 
590= 
600= 
610= 
620= 
630= 
640= 
650= 


PROGRAM  D I  SK£4  <  I  NpljT .  OUTPUT  »  TAPE5=  I  NPiJT .  ?APE6*0UTPUT> 

C  FAR  FIELD  SCATTERING  BY  A  CIRCULAR  METALLIC  DISK 

D I  MENS  I  ON  LABEL  <5)  .  VAR  <  5) 

COMPLEX  F4 .  PS  I  .  W .  VO  -  YET  AO  •  U  .  X  .Z.ZA»2B>  ZC >  ZD 
1 .  I X .  EF'ART «  EPftRP *  EF'ERT .  EF'ERP  •  ESNT .  ESNP 
COMMON  EI6  <5 0>  >  D  <50 . 5 0>  *  DNEG  <5 0>  >  P 1  <5 0>  >  F4  <5 0> 

1 » SO  <50>  . P  <50>  *  SETA  <50>  » SETAO  <50>  > PS  I  <50>  » W  <50>  *  YO  <5G> 

1 *  YET  AO  <5  0>  > IJ  <5  0>  » X  <5  0>  » CMPH I  <5  0>  »  SMPh  I  <5  0) 

DATA  LABEL  <1>  ■'1  OH  KA 
DATA  LABEL  <2)  s 1  OH  THETA  I  ✓ 

DATA  LABEL  <3)  1  OH  F'DL  ✓ 

DATA  LABEL  <4>,'l  OH  THETA  S  ✓ 

DATA  LABEL  <5>  ■"'1  OH  PHI  S  / 

IX=<0. »  1.) 

WRITE  <6.  11) 

11  FORMAT < 1 X. /// . IX. "SCATTERING  BY  A  METALLIC  CIRCULAR  DISK 

1 5X."  CHDDGE  —  VERSION  12/17/78) ") 

WR I TE  <6  .  £6) 

£6  FORMAT  <  1 X »  /// ,  1 X  .  "  <T YPE  KA=  0  TO  STOP.  PROGRAM)  “  .  /  > 

1 IX. " (TYPE  KA=- 1  FDR  A  DESPIPTION  OF  THE". 

1 "  PARAMETERS) " . /. IX. " (NORMALIZATION:  ESC AT= A*E I NC^ENOPM 
1  "  /  <£*P>  ♦EXP  <- J^>P  :•  )  " .  / ,  1 X »  "  CALL  ANGLES  I N  DEGREES)  ”  ) 

4  NOC=l 

INDEX=1 
WRITE <6» £7) 

£7  FORMAT  <1X. ///? IX. " 1 .  KA” . 14X. ”=  ”) 

READ  <5. ♦)  VAR  < 1 ) 

IFCVARCl)  .EQ.-DGO  TO  41 
I F  CVAR <1 ) . LE . 0 . ) GO  TO  5 
WRITE <6. 28) 

£8  FORMAT  <1X» ”£.  THETA  INCIDENT  =  ”) 

READ  <5 . ♦)  VAR <£) 

WRITE <6. 29) 

£9  FORMAT C1X. "3.  POLARIZATION  =  ”) 

READ  <5. ♦)  VAR (3) 

WR  I  TE  <6 .  3  0) 

30  FORMAT <1X. "4.  THETA  SCATTERED  =  “> 

READ  <5.  ♦)  VARC4) 

WRITE <6. 31) 

31  FORMAT C1X. "5.  PHI  SCATTERED  =  ”) 

READ  <5.  ♦)  VARC5) 

WRITE  <6. 3) 

3  FORMAT ax.yy. IX. "WHICH  VARIABLE  IS  TO  BE  INCREMENTED?") 
READ  <5. ♦)  NVAR 

I F  < CNVAR . LE . 0) . OR . CNVAP . GT . 5) ) GO  TO  £3 
WRITE <6. 21) 

£1  FORMAT <1X. "TYPE  NUMBER  DF  CASES:") 

READ  <5. ♦)  NOC 
WRITE <6. 32) 

32  FORMAT <1X. "WHAT  IS  THE  INCREMENT?") 

READ  <5 » ♦)  VINCRE 

£3  CONTINUE 

WR I TE  <6 . £4)  LABEL  CNVAR) 

24  FORMAT  < 1 X . // . A 1 0 . 4X . " CROSS  SECT I ON " . ££X . " E  NORM " .  / . 

1 1 3X .  " S I GMRy  CP I ♦A«*2) " . 1 1 X . " THETA " . 1 5X . "PH I" . / . 1 3X . 

1 "THETA" . 7X> "PHI ” * SX. "MAG".5X. "PHASE". 6X. "MAG".5X. 
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66  0= 

1" PHASE"*/} 

6?  0= 

6 

C=VAP  C  1  > 

680= 

MMMRX=45 

68  0= 

MMMAX=45 

70  0= 

IRRMAX=45 

710= 

THEO=VAP  ■•£>  ♦S .  1 4 1 59/  1 8  0 

720= 

ET80=C0S  CTHEO) 

730= 

THE=VAP  ■:4>  *3.  14159/180 

740= 

ETR=COS  CTHE} 

750= 

PH  I  =  VHP  05>  ♦3.141 59/ 1 8  0 

760= 

RLF=VRR  <3> *3 .141 59/ 1 8 0 

770= 

CRLF=COS CPLF} 

78  0= 

SRLF=SIM  CALF} 

790= 

34 

DO  10  MM= 1 * MMMRX 

800= 

M=MM— 1 

810= 

10=0 

82  0= 

CALL  SMNO CM - MMMRX} 

830= 

CRLL  DEE I GN <C *  M * MMMRX} 

84  0= 

CALL  OBCOFM <C » M * MMMRX , MMMRX » I RPMRX>  . 

850= 

CRLL  DMEGM  CC * M> MMMRX} 

860= 

CRLL  OERRD  CC *  M «  I Q *  MMMRX *  MMMRX  -  I RRMRX} 

870= 

IF  CIO. EC. 1>  GO  TO  10 

880= 

CRLL  FPSI CM* MMMRX} 

890= 

CRLL  POLYM  C ETfiO *  M * I RRMRX} 

90  0= 

CRLL  0BRN6CNNMAX* I RRMRX} 

910= 

DO  7  1=1 > MMMRX 

920= 

7 

SETRO  CI>  =SETR  CD 

930= 

CRLL  FM CM » MMMRX} 

940= 

CALL  POLYM  CETfi * M * I RRMRX} 

950= 

CALL  0BRM6 CNMMRX » I RRMRX} 

96  0= 

CRLL  FY CM* MMMRX} 

970= 

10 

COMTIMUE 

980= 

DO  13  MM=1 * MMMRX 

990= 

CRLL  FXIJ  CMM  *  MMMRX} 

1000= 

IF CMMMRX.LT. MM} 60  TO  13 

1010= 

CALL  CSPHI CMM* PHI} 

1020= 

13 

COMTIMUE 

1030= 

CRLL  FZ  CMMMAX*  Z*  ZR*  ZB*  ZC*  ZD} 

104  0= 

EPART=ETA* •-2^Z^CMPHI C2} +ZR} 

1050= 

EPRRP=-£^Z^SMPHI C£}+ZB 

1060= 

EPEPT=ETR^  C£^Z^SMPHI  •:'£}  -ZC} 

1070= 

EPEPP=-2^Z^CMPHI C£} +ZD 

1080= 

IF CETRO. EO. 0. } GO  TO  2 

1090= 

ESMT=2MX^  CCRLF^EPRPT/ETRO+SRLF^EPERT} 

1100= 

ESMP=2MX^  CCRLF^EPRRP/ETRO+SRLF^EPERP} 

1110= 

GO  TO  8 

1120= 

2 

ESMT =£♦  IX^SRLF^EF’ERT/C 

1130= 

E$NP=2MX*SRLF#EPERP/C 

1140= 

8 

EMRGT=CRBS ’ESMT} 

1150= 

EMRGP=CRBS CESMP} 

1160= 

S I 6THE=EMRGT ♦EMRGT 

1170= 

S I GPH I =EMRGP*EMRGP 

1180= 

E1=REAL  CESMT} 

1190= 

E2=R I MRG  CESMT} 

1200= 

IFCE1.EO.  0. } GO  TO  16 
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1210= 

1220= 

1230= 

1240= 

1250= 

126  0= 

16 

1270= 

1280= 

17 

1290= 

1300= 

1310= 

1320= 

1330= 

134  0= 

1350= 

1360= 

18 

1370= 

1 38  0= 

19 

1390= 

1400= 

1410= 

£5 

1420= 

37 

1 43  0= 

144  0= 

1450= 

1460= 

5 

1470= 

41 

1480= 

40 

"?/? 

1490= 

?1X? 

1500= 

1510= 

1520= 

1530= 

154  0= 
1550= 
1560= 
1570= 
1580= 
1590=  42 
1600= 
1610= 

1 62  0=C 

1 63  0=C 
1640=C 
1 65  0=C 
1660=C 
1670= 
1680= 
1690= 
1700= 
1710= 
1720= 
1730= 
1740=  2 
1750=  1 


APG=E£-E1 

EF'HAT  =  1 8  0-3 .141  59^ATAN  <  ARG 

IF  <  <E1 .  LT.  0.  >  .  Fhlt.  <E£. GT.  0.  >  .>  EF'HAT =EF'HAT +180 
I F  <  <E1.LT.  0 .  >  .  FiND .  <E£ .  LT .  0.  >  >  EF’HAT=EPHAT- 1 8  0 
60  TO  17 
EPHAT=90 

IFCE2.LT.  0. >EPHAT=-90 
E1=REAL<ESNP> 

E2=AIMAGCESNP> 

IF  <E1 .  EG1.  0.  >  GO  TO  18 
ARG=E2'E1 

EPHAF- 1 8  0--3 .141 59*ATAN  <ARG> 

IF  <  <E1 .  wT.  0.  >  .  FIND .  <E£ .  GT.  0.  >  >  EF'HRP=EPHAP+ 1 8 0 
IF  <  <E1 .  LT.  0.  >  .  RND.  CEE.  LT.  0.  >  >  EPHRF-EF'HAF-180  1 

GO  TO  19  j 

EPHRP=90  I 

IFCE2.LT.  0.:>EPHRF=-90 

I F  CNVRP . EQ . 0> NVRP= 1  j 

MR I TE  <6  ? £5>  VRR CNVRR> > $ IGTHE • 3 1 GPH I ? EMAGT  *  EF'HAT  ? 
1EMAGP?EPHAP 

FORMAT  <  1 X  ?  F7 .  £  ?  2  <  1 X  ?  E 1  0 .  3>  ?  £  <  1 X  ?  E 1  0 .  3  ?  1X»F7.2>> 
IFCINDEX.EQ.NOOGD  TO  4  j 

I NDEX= I NDEX+ 1  i 

VRP <NVAR>  =VRP <NVAR> +V I NCRE  1 

GO  TO  6 
CALL  EXIT 
WRITE  <6*  4  0> 

FORMAT < IX? IX? "THE  DISK  OF  RADIUS  A  LIES  IN  THE  X-Y  PLANE 

1 1 X  ?  " CENTERED  AT  THE  ORIGIN.  THE  CONVENTIONAL  <R?  THETA?  "  ?  / 

1 " PH I >  COORDINATE  SYSTEM  IS  USED  IN  THE  FAR  FIELD.  ">/»lX. 

1  "THE  PLANE  OF  INCIDENCE  OF  THE  PLANE  WAVE  IS  THE" » ■-'?  IX? 
1"X-Z  CPHI  =  0>  PLANE.  THE  POLARIZATION  ANGLE  <POLV‘ ? /? IX? 

1  "OF  EINC  IS  MEASURED  FROM  THE  PLANE  OF  INCIDENCE  "?■-'?  IX? 

1 "  IN  THE  PHI-DIRECTION.  I.E.?  PDL=0  IS  THE  PARALLEL"  ?  ✓"?  IX? 

1 " CTHETA)  CASE  AND  POL=90  IS  THE  PERPENDICULAR  CPHI>"?/?1X? 
1"CASE. "?£X? "THE  RESULT  IS  A  SOLUTION  OF  THE  RIGOROUS" ? s* 

1 1 X ?  " E I GENFIJNCT I  ON  SCATTER  I NG  PROBLEM .  "  > 

GO  TO  4 
CONTINUE 
GO  TO  4 
END 

OBLATE  SPHEROIDAL  ANGULAR  FUNCTION? S? OF  ARGUMENT  0? 

ORDERM?  N5  WITH  N-M  EVEN? UP  TO  ORDER  N=M+£^NNMAX-£. 

CEQUAL  TO  PMN C  0> > 

SUBPOUT 1 NE  SMNO CM  ?  NNMAX> 

COMPLEX  F4 

COMMON  EIG  C50>  ?  D  <50?  50>  ?  DNEG  C50>  ?  R1  <50>  ?  F4  <50> 

1  ?  SO  C50>  i 

S0C1>=1 

IF  CM. EQ. 0. >  GD  TO  1 

DO  2  MM= 1 ?  M 

SO  <  1  >  =  C2^MM- 1 *SO  C 1  > 

DO  3  NN=1 ?  NNMAX 
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1760=  n=£+CNN-1>+M 

1 77  0=  3  S  O  CNN+ 1 > a- CN+M+ 1 > +SO CNN> ✓ CN-M+£> 

1780=  RETURN 

1790=  END 

180  0=C 

181 0=C  SIN  FIND  COSINE  FUNCTIONS  OF  PHI 

1 88  0=C 

1 83  0=  SUBROUT I NE  CSPH I  CNN . PH I > 

184  0=  COMPLEX  F4 . PS I  *  W ,  YO  *  YETfiO » U  *  X 

1 85  0=  COMMON  E 1 6  >-5 0>  - 1>  C5 0 » 5  0>  .  DNEG  <  5  0>  »  R 1  C5  0>  »  F4  <5  0> 

1 86  0=  1  r  SO  C5  0>  » P  C5  0>  > SETH  <5  0>  » SETfiO  C50>  •  PS  I  C5  0>  t  W  <5 0)  , YO  C5  0) 

1 87  0=  1 »  YETfiO  <5  0>  »  U  <5  0>  » X  C5  0>  » CMPH I  C5  0>  »  SMPH I  <5  0> 

1880=  M=MM-1 

1 89  0=  CMPH I CMM> =COS CM+PH I > 

1900=  SMPH I CMM> =S I N  CM+PH I > 

1910=  RETURN 

1980=  END 

1 93  0=C 

1 94 0=C  OBLRTE  SPHEROIDAL  EIGENVALUES  OF  ARGUMENT  C»  ORDER  M»N 

1 95 0=C  WITH  N-M  EVEN  UP  TO  ORDER  N=M+£+NNMfiX-£ 

1 96  0=C 

1970=  SUBROUTINE  OBE I GN <C >  M  > NNMflX> 

1980=  COMMON  EI6C50) 

1 99  0=  D I  MENS  I  ON  I P  C5  0>  *  ALPHA  <5  0>  »  BETA  C5  0>  *  P  C5  0> 

3000=  4  CONTINUE  , 

2010=  M£=2+M 

2020=  C£=C+C 

2030=  ACC=1 • 0E-05 

£040=  NN2=NNMAX+£ 

2050=  NI=NN£+1 

2060=  P  C1>=1 

2070=  IPC1>  =  1 

2080=  DD  £  I  QG'=  1  ?  NN2 

2090=  IV=£+I00-1 

2100=  IW=M£+£+IGO 

2110=  I X=M2+4+I QQ— 1 

2120=  ALPHA  C I  QQ>  =  CC2+  CM2+  <£♦  I V- 1  >  +£♦  I V+  <  I V- 1  >  - 1  >  >  ✓  C I X+  C I  X-4)  > 

2130=  1 - CM+ I V- 1 > ♦ Cl V+M> 

2140=  £  BETA  C I QO+ 1  >  =C£/ 1  X+SORT  C I V+  C I V+ 1  >  ♦  I W+  C I W- 1  >  ✓  C I  X+  I X-4 .  0>  > 

2150=  BETA  CNN2+ 1 >  =  0 . 

£160=  BO=ABS  CALF'HA  C 1  >  >  +ABS  CBETA  C£>  > 

2170=  DO  3  I OQ=£ » NN2 

£180=  flO=ABS '  BETA  CIQO) >  +AES CALPHfi CIQO> ) +ABS CBETA Cl 00+ 1> > 

2 1 9 0=  BETA  C 1 00> =BETA C I QCO  +BETA C I QO> 

£200=  IF CAO. GT. BO>  BO=AO 

2210=  3  CONTINUE 

£220=  AO=-BD 

2230=  BOI=BO 

£240=  13  CONTINUE 

2250=  BO=BOI 

2260=  DO  £0  100=1 j  NNMfiX 

££70=  N=£+I0Q-2+M 

£280=  fi=A0 

2290=  B=BO 

£300=  IERR=-1 


£1 


£31  0= 
£320= 
£33  0= 
£34  0= 
£350= 

c)  t*  0= 

£370= 
£330= 
£39  0= 
£400= 
£410= 
£4£0= 
£430= 
£440= 
£450= 
2460= 
2470= 
2480= 
£490= 
2500= 
£510= 
£520= 
£530= 
£540= 
£550= 
£560= 
£570= 
£580= 
£590= 
£600= 
£610= 
2620= 
£630= 
£64  0= 
2650= 
266  0= 
£670= 
2680= 
2690= 
2700= 
£710= 
£720= 
2730= 
£74  0= 
2750= 
£760= 
2770= 
£780= 
2790= 
£800= 
£810= 
£820= 
2830= 
2840= 
£850= 


50 


41 

4£ 

40 


5 

12 


C' 

14 

q 

10 

11 


15 


16 


£4 
7  rtf i 
£0 


C 

C 

C 

C 

C 


4 


113=0 

CO=  <R+B>  ■■•"£ 

I F  ('CD  >  5  0 *  ££  *  5  0 
ERR=  >:F-A>  .'-PB3  <CO> 

IEPP=IEPP+1 

IF  a  EPF-6 n:>  4 O'.  41 -.41 

l,IF:ITE<6«4£>  N 

FORMAT  <£X»  "ITERATIONS  EXCEEDED  FDR  EIGENVALUE  "*I3> 

GD  TO  700 

I F  <EPP-ACC>  £4 *  £4  •  ££ 

P  <£>  =RLF‘HA  <  1  >  —CD 
DO  5  1=3* N I 

P  <1  >  =  < ALPHA  <1  - 1  >  -CO-BETA  <1  - 1  >  ♦  <!P  <1  -£>  /P  <1  - 1  >  >  >  +P  <1  - 1  > 
PMAG=AE3  <P  <I> > 

IF  <PMAG. GT. 1 . OE+33)  GO  TD  7 

CONTINUE 

CONTINUE 

DO  6  I =£ * N I 

IF<P<I>>  1 4 !>  8 j  9 

IF  <P  <I-1>  >9*9* 14 

I P  < I >  =— 1 

GO  TO  10 

ipci>=i 

IF<IP<I>-IP<I-1>>6*  11*6 

113=113+1 

CONTINUE 

I F  < 1 1 S- 1 QQ> 16* 15* 15 

R=CO 

GO  TO  21 

B=CO 

GO  TO  £1 

BO=CO 

EIG  <IQ0>=-CO 
CONTINUE 
RETURN 
NNMAX= I —4 
N I =NNMRX+3 
GO  TO  4 
END 

OBLATE  SPHEROIDAL  EIGENFUNCTION  EXPANSION  COEFFICIENTS 
OF  ARGUMENT  C*  ORDER  M *  N *  R *  WITH  N-M  EVEN*  UP  TO  ORDER 
N=M+£+NNMAX-£  >  R=£+ 1 RRMAX+2 

SUBPOUT I NE  DBCOFN < C *  M  *  MMMAX *  NNMRX *  I RPMAX> 

COMMON  E I G  <5 0>  *  D  <5  0  *  5  0> 

DIMENSION  DP <50> 

C£=C+C 

MM=M+1 

DO  1  NN=1 *  NNMAX 
N=M+£+NN-£ 

DP  < I RRM AX+3>  =  0 
BP  <IRRMAX+£>  =  1 . 0E-30 
D  CNN* 1>  =  0 
D  <NN*  £>  =  1 
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£860= 

J  J=  ■  N-M>  •••'£+ 1 

£87  0= 

DO  107  LL=1 * IRRMRX 

£880= 

L=LL- 1 

£890= 

I F  CLL . GE . J -0  L= I PFMRX+ J J-LL 

£900= 

IF'=£^L 

£910= 

IPM=M+IR 

£9£0= 

RP=  <M+ 1  PM+£>  ♦  <  M+ 1 PM+ 1  >  ♦C£^  C  <£♦  I  RM+3>  ♦  <£♦  I  RM+5>  > 

£930= 

P.p=  •:'£♦  I PM^  C  I PM+1  >  -£»I1»H- 1  >  ♦C£^  C  <£♦  I  PM- 1  >  ♦ 

£94  0= 

1  <£♦  I  PM+3>  >  - 1  PM*  C  I PM+ 1  > 

£•=>=:  o= 

CP=  I R>  <I P- 1  >  ♦•:£-••  <  <£♦  I PM-3>  ♦  <£♦  I PM- 1  > > 

£960= 

IF  <LL-JJ> 1 05*  1 06 *  1 06 

£970= 

105 

D  CNN *  L+3> =- CCR*D CNN *  L+l>  +  CBR+E 1 6  CNN> > ♦D  CNN * L+£>  > /fiR 

£98  0= 

DMRG=RBS CD CNN * L+3> > 

29i,=i  o= 

IFCEMRG.GT.  1. 0E+30>  GO  TO  3 

3000= 

GO  TO  107 

3010= 

1  06 

DP CL+1> =- CRP^DP CL+3>  +  CBR+E I G  <NN> > ♦DP CL+£>  >  xCR 

30£0= 

DMRG=RBS<DPCL+1>> 

3030= 

IF  CDMRG. GT. 1 . 0E+30>  GO  TO  3 

3040= 

107 

CONTINUE 

3050= 

DL=RBS CD  CNN *  JJ+ 1 > > 

3060= 

DL=RL0G1 0 CEL) 

3070= 

DLP=RES CDPCJJ+1)) 

3080= 

DLF'=RLDG1 0  CDLP) 

3090= 

DL=RES CDL) 

3100= 

DLP=RPS CDLP) 

3110= 

BL=DL+DLP 

31£0= 

IF  CDL. GT. 30.  >  GO  TO  5 

3130= 

CON=D CNN  j  J J+ 1 > ✓DP CJJ+1 > 

314  0= 

RCON=flBS CCON> 

3150= 

IF CfiCON. GE. 1 . 0E+3£)  GO  TD  £ 

3160= 

DO  118  J=JJ* IPRMRX 

3170= 

118 

D  CNN*  .J+£>  =CON+DP  C.J+£> 

3180= 

F=  1 

3190= 

IF  CM) 198* 198* 199 

3£  0  0= 

199 

DO  110  1=1 *M 

3£  1  0= 

110 

f=f+cm+i:> 

3££  0= 

198 

SUM=  0 

3£3  0= 

MMX= I RRMRX+ 1 

3£4  0= 

DO  113  1=1* MMX 

3£5  0= 

IR=£+I 

3£6  0= 

SUM=SUM+F+D  CNN  *  I  + 1  > 

3£70= 

IF  CI-JJ)  113*197*113 

3£80= 

197 

FNM=F 

3£90= 

113 

F= C-F+ C I R+£+M- 1 )  )/IR 

3300= 

RLF=FNM--'SUM 

3310= 

DO  114  1=1* MMX 

33£0= 

D  CNN  *  I )  =RLF+D CNN  *  I  + 1 > 

3330= 

114 

CONTINUE 

3340= 

1 

CONTINUE 

3350= 

RETURN 

336  0= 

3 

IRPMRX=LL-1 

3370= 

GO  TO  4 

3380= 

£ 

I RRMRX= I RPMRX- 1 

3390= 

GO  T04 

3400= 

5 

NNMRX=NN-1 
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3410=  RETURN 

3420=  ENI' 

343 0=C 

344 0-C  NEGATIVE  D  COEFFICIENT  SUBROUTINE 

345 0=C 

346  0=  SUBPOUT I NE  DNE6N < C - M * NNMRX> 

3470=  COMMON  E I G  <5  0>  *  D  <5  0 >  5  0>  *  DNEG  <50> 

3480=  DO  4  NN= 1 > NNMRX 

3490=  C£=C^C 

3500=  I F  <M . GE . 1 >  GO  TD  £ 

3510=  HO  5  11  =  1*  NNMAX 

35£ 0=  5  DNEG < 1 1 > »D  < 1 1  *  1 > 

3530=  GO  TO  3 

3540=  £  B 1=1.0 
3550=  B£=0. 0 

3560=  EI=EIGCNN> 

3570=  HD  1  IPP=1*M 

3580=  I  R=£'4 1  PR-£^M-£ 

359  0=  RR=  <£-*M+ 1  R+£>  ♦  •:£^M+ 1 R+ 1  >  +C£s  <  <£4M+£4 1  R+3> 

36  00=  1 ♦ <£^M+£*I R+5> > 

36 1 0=  BR= CM+ I R> ♦ <M+ I P+ 1 > -E I - <£♦ <M+ I P> ♦ <M+ I R+ 1 > -£*M*M- 1 > 

36£  0=  1  ♦C£.-"  <  <£4M+£4 1 R- 1  >  ♦  <2H1+£+ 1 R+3)  > 

363  0=  CR=  <  I  R>  ♦  <  I R- 1  >  ♦C£/  <  <£♦«+£♦  I  R-3>  ♦  C£+t1+£+ 1 R- 1 )  ) 

3640=  B3=B£ 

3650=  B£=B 1 

366  0=  1  B 1  =  <BP4B£-CR4B3>  'AR 

3670=  R=B  <NN» 1) /B1 

3680=  DNEG <NN> =fl 

369  0=  DDM=ABS  <DNEG  <NN>  > 

3700=  IF <DDM. LT. 1 . 0E-35>  GO  TO  6 

3710=  4  CONTINUE 

37£  0=  3  RETURN 

3730=  6  NNMAX=NN-1 

3740=  RETURN 

3750=  END 

376 0=C 

377 0=C  OELRTE  SPHEROIDAL  RADIAL  FUNCTION  R  <4>  OF  ARGUMENT  C 
378 0=C  ORDER  M * N 5  WITH  N-M  EVEN*  UP  TO  ORDER  N=M+£*NNMRX-£. 
379 0=C  ALSO  NORMALIZATION  FUNCTION*  N. 

3800=C 

38 1 0=  SUBROUT I NE  OBR AD  <C*M* IQ* MMMAX >  NNMAX *  I RRMRX> 

38£  0=  COMMON  E 1 6  <5  0>  *  D  <5 0  *  5  0>  *  DNEG  <5 0>  >  R 1  <5 0)  *  F4  <5  0> 

3830=  COMPLEX  IX*R4*F4 

3840=  IX=<0.0*1.0> 

3850=  EEAC=1 

3860=  EAC=1 

3870=  EAC£=1 

3880=  GRO=l 

3890=  IF  <M.  EQ.  0>  GO  TO  £0 

3900=  MRXM=M+ 1 

3910=  DO  19  MM=£  *  MAXM 

39£0=  IM=MM-1 

3930=  GRD=  <£♦ I M- 1 > ♦ ££♦ I M> ♦GPO 

3940=  EEAC=  <£♦ I M- 1 > ♦EEAC 

395  0=  E AC£=  <£♦ I M- 1 > ♦ <£♦ I M> ♦EACE 
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3960=  19 

ERC=IM«€RC 

3970= 

I F  >  EEhC . GE .  1 . 0E+1 7>  GO  T 

393  0= 

IF  '  E H L c‘  •  IrE  .  1 .  uE+30>  GO  T 

3990=  £0 

DO  17  NN= 1 ?  NNMRK 

4  000= 

N=£^CNN-1>+M 

401  0= 

SUM*  0 

4  02  0= 

bP=GRO 

4030= 

EN0PM=0. 

404  0= 

DO  IS  NR=1 »  IPRMfiX 

4  050= 

IR=£^ CNR-1> 

4  06  0= 

SUMP=GR^D  CNN*  NR> 

4070= 

S UN = SUM + SUMP 

4030= 

DMRG=RBS  CD  CNN  *  NR>  > 

4090= 

I F  C'DMRG .  LT .  1 . 0E-3 0>  GO  T 

4100= 

E N D R M P = £ ♦  G R ♦  CD  CNN*  NR>  > 

4110=  £ 

E  N  0  P  M = E  N  O  R  M + E  N  0  P  M  P 

4 1  £  0= 

GRMRG=RBS  C  G R > 

4130= 

IF  CGRMRG. GT. 1 . 0E+30>  GD 

4140= 

6R=  C  I  R+£-*M+ 1 .  >  /  C  I R+ 1 .  >  ♦  < 

4150=  18 

CONTINUE 

4160= 

R 1  <NN>  =  C- 1  >  ♦♦  CNN- 1  >  ♦£♦^1' 

4170= 

R£=  C-l  >  ♦♦  CNN- 1  >  ♦  C£*M—  1  >  ♦( 

4130= 

1  ♦  CE  E  R  C  ♦♦£.-'  S  U  M  >  ♦£**M.-DNEG 

4190= 

PR£=PBS CR£> 

4200* 

IF  CRR£. GE. 1 . 0E+30>  GD  TO 

4£1  0= 

EERC=  CN+M+ 1  >  ♦EEPC  --  <N-M+£‘; 

4££0= 

I F  CEERC . GE . 1 . 0E+ 1 7>  NNMR! 

4£30= 

R4=R1  CNN>-I»R£ 

4£40= 

R5=RBS <ENORM> 

4250* 

R 1 =RLOG 1 0  <R5> 

4£6  0= 

R4=RBS <R2> 

4£7  0= 

Rc-RLOGl 0  <R4> 

4£8  0= 

R3=R 1 +R£ 

4£9  0= 

IFCfi3.GT.30.>  GO  TO  3 

4300= 

R5=RBSCR1 <NN>> 

4310= 

Rl=RLOGl 0 <R5> 

43£0= 

R3=RES CR1> +RBS CR£> 

4330= 

IFCfl3.GT.30.>  GO  TD  3 

4340= 

F4 CNN> = 1 / CEN0RM^R4> 

4350=  17 

CONTINUE 

4360= 

RETURN 

4370=  1 

ENORMP* 0 

4330= 

GO  TO  £ 

439  0*  3 

F4  CNN>* CO. *  0.> 

4400= 

GO  TO  17 

4410=  5 

GR=0 

44£  0= 

GO  TO  IS 

4430=  4 

MMNfiX=M 

4440= 

IQ=1 

4450= 

RETURN 

4460= 

END 

447 0=C 

443 0=C 

RSI  FUNCTION 

449 0=C 

4500= 

SUBROUTINE  FPSI CM* NNMRX> 

4 

4 


I  R+£^M+1) 


TD 


cnn> 
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45 1 0=  COMPLEX  I X . F4 , RS I 

45c'  0=  COMMON  E I G  '  5  0>  » D  <  5  0 «  5  CO  .  DNEG  <5  CO  •  P 1  <5  CO  >  F4  C5  0> 

453  0=  1 » SO  <5  CO  » P  <5  0)  *  SETA  <5  CO  *  SET  AO  <5  CO  *  PS  I  <5  0) 

4540=  IX=C0.»1.> 

4550=  MM=M+ 1 

456  0=  PS  I  CMM>  =  C  0 .  »  0 .  > 

4570=  DO  1  NN= 1 >  NNMAX 

4580=  N=M+£^NN-2 

4580=  F'S  I  CMM>  =PS I  CMM>  +IX^  <N>  ♦SO  CNN>  ♦♦c>F4  CNN> 

4600=  1  CONTINUE 
4610=  RETURN 

4620=  END 

463 0=C 

464 0=C  ASSOCIATED  LEGENDRE  POL YDM I ALS >  OF' >  DF  ARGUMENT  ETAD? 
465 0=C  ORDERPMjN;  WITH  N-M  EVEN?  DF  TO  ORDER  N=M+£^NNMAX-£. 

466 0=C 

467  0=  SUBROUT I NE  POL VN  CETAO » M • NNMAX> 

4680=  DIMENSION  PPC3> 

4690=  COMPLEX  F4 

47 0  0=  COMMON  E I G  <5  CO  «  D  <  50  5  0>  DNEG  <5  CO  *  R 1  <5  CO  ->  F4  <5  0> 

4710=  1  >  SO  <5  0>  » P  <5  CO 

4720=  SO=SQRT Cl-ETAO^ETAO> 

4730=  F‘P  <  1  >  =  0 

474  0=  PP  <£>  =1 

4750=  IF  CM.  EG!.  0>  GO  TO  1 

4760=  DO  2  L=  1  *  M 

4770=  2  PP  <£>  =  C2^L-1>  ♦SCf^PP  <2> 

4780=  1  F'C1>=PPC2> 

4790=  DO  3  NN=£»  NNMAX 

4800=  N=M+2^NN-3 

4810=  DO  4  L= 1 ?  £ 

482  0=  PP  <3>  =  < C2^N- 1 > ♦ETAO^PP C2> - CN+M- 1 > ^PP  C 1 >  > ✓ CN-M> 

4830=  N=N+ 1 

484  0=  PPC1)=PP<2> 

4850=  4  PPC£>=PP<3> 

4860=  3  P CNN> =PP <3> 

4870=  RETURN 

4880=  END 

489 0=C 

4900=C  OBLATE  SPHEROIDAL  ANGULAR  FUNCTIONS-.  S»  OF  ARGUMENTS 

491 0=C  C  AND  ETAO;  ORDER  M»N5  WITH  N-M  EVEN;  UP  TO  ORDER 

492 0=C  N=M+£^NNMAX-£. 

493 0=C 

494  0=  SUBROUT I NE  OB ANG CNNMAX » I PRMAX> 

4950=  COMPLEX  F4 

496  0=  COMMON  E I G  <5  CO  »  D  <5  0 »  5  0>  »  DNEG  C5  0>  » R 1  <5 0>  »  F4  C5  CO 

497 0=  1 » SO  <  5  CO  ,  P  <5  CO  >  SETA  <5  0>  *  SETAO  C50> 

4980=  DO  1  NN= 1  *  NNMAX 

4990=  SETA  <NN> =0 

5000=  DO  £  IRR=1  ->  IRRMAX 

5010=  £  SETA  <NN> =SETA CNN>  +D CNN » I RR> ♦P  C I RR> 

5020=  1  CONTINUE 

5030=  RETURN 

5040=  END 

5050=C 


82 


5060=0  W  FUNCTION 

5070=0 

5060=  SUFROUTINE  FW  (M*  NNMfiX'1 

5090=  COMPLEX  IX*  F4,  l.i.  PCI 

51  00=  COMMON  E  16  ■  50 >  *  D  1  50*  5 0>  *  DNEG  (50‘>  *  F'l  <50>  *  F4  <'50> 

5110=  1  *  SO  1 -'5 O')  *  p  <50>  •  SETR  <50  '*  *  SETRO  (  5 0 *  PS  I  (  5  0  >  *  W  (50) 

5160=  IX= >: 0. , 1.) 

5130=  MM=M+ 1 

514  0=  1.1 (MM>  =  ( 0 . *  0 . > 

5150=  DO  1  NN= 1  *  NNMfiX 

5160=  N=M+£^NN-£ 

5170=  1  t.l  (MM>  =l.l  <MM>  + 1  X^N^SETfiD  (NfO  ♦SO  <NN>  ^F4  <NN> 

5180=  RETURN 

5190=  END 

5£  0  0=C 

5£ 1 0=C  Y  FUNCTION 

5££ 0=C 


523  0= 

SUBROUT I NE  FY  (M  *  NNMfiX) 

5£40= 

COMPLEX  F4 *  YO *  YETfiO *  PS I  *  W 

525  0= 

COMMON  E I G  (5  0>  *  D  (5 0  *  5  0>  *  DNEG  (5  0>  *  R 1  (5  0>  * 

F4(50> 

5260= 

1  *  SO  (5  0>  *  P  (5 0>  *  SETR  (5 0)  *  SETRO  (5  0>  *  PS  I  (5  0!> 

,  I.J  (5  CO  *  YO  (50) 

5270= 

1*  YETfiO  (50) 

528  0= 

MM=M+1 

5290= 

YO  (MM>  =  (0.  *  0.  > 

5300= 

YETfiO  (MM>  =  (0.  *  0.> 

5310= 

DO  1  NN=1* NNMfiX 

5320= 

N=M+£+NN-£ 

533  0= 

YO  (MM> =YO  (MM>  +  (- 1 > ♦♦N+P 1 (NN> ♦SO (NN> ♦SETR 

(NN> +F4 (NN> 

5340=  1 

YETfiO  (MM>  =YETRO  (MM>  +  <- 1  >  ♦♦((♦R 1  (NN>  ♦SETRO 

(NN>  ♦SETfi(NN)  ^F4  (Nf 

5350=  RETURN 

5360=  END 

537 0=C 

538 0=C  S  RND  U  FUNCTIONS 

539 0=C 

54  0  0=  SUBROUT I NE  FXU (MM *  MMMfiX> 

54 1 0=  COMPLEX  I X *F4*PSI» W *  YD  *  YETfiO *  U *  X *  PT 

54£  0=  COMMON  E 1 6  (5  0>  *  D  (5  0  *  5  0>  *  DM  EG  (5  0>  *  R 1  <5  0>  *  F4  <5  0) 

543  0=  1  *  SO  (5  0>  .*  P  (5  0>  >  SETfi  <5  0>  >  SETRD  <5  0)  *  PS  I  <5  0>  *  W  <5  0>  *  YO  <5  0> 

544  o=  1 ,  YETfiO  <5 0>  *  U  (5  0>  *  X  <5  0> 

5450=  IX=(0.*1.> 

5460=  IF  (MM. EQ. 1>  GO  TO  1 

5470=  I F  (MM . EO . MMMRX)  GO  TO  £ 

548 0=  PT=PS I (MM- 1 > ♦PS I (MM+ 1 > 

5490=  PTT=CfiBS (PT> 

5500=  IF (PTT. EQ. 0. >  GO  TO  3 

55 1  0=  U  (MM>  =£♦  I X^ <MM— £>  ♦  (l,.i  (MM- 1  >  +1.1  (MM+ 1>> ✓  (PS I  (MM- 1  >  +PS I  (MM+ 1 )  > 

55£ 0=  X  (MM>  =£♦  I  Y.++  CMM-2)  ♦  (W  (MM- 1  >  -W  (MM+ 1  >  .)  ✓  (PS  I  (MM- 1  >  +PS I  (MM+ 1  >  • 

5530=  RETURN 

554  0=  1  U  ( 1  >  =- 1  X^W  (£>  z'PS  I  (£> 

5550=  X  ( 1 .)  =  ( 0 .  *  0 .  > 

5560=  RETURN 

557  0=  £  U  (MM>  =  I X++  (MM-£>  ♦£♦1.1  (MM- 1  >  /PS I  (MM- 1 ) 

5580=  X  (MM) =U  (MM) 

5590=  RETURN 

5600=  3  MMMfiX=MM-l 
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5610=  RETURN 

5630=  END 

56:-:  n=r 

564 0=C  Z  FUNCTIONS 

565  0=C 

566  0=  SUBPOUT  I NE  FZ  CMMMRX  •  Z  *  ZR ?  ZB ?  ZC  ?  ZB> 

5670=  COMPLEX  IX?  10?  Z?  ZR  ?  ZB?  ZC?  ZD?  F4 »  F'SI  ?  W?  VO?  YET  RO  ?  U?  X 

568  0=  COMMON  E I G  <  5  0>  *  D  <5  0  ?  5  0>  ?  DNEG  5  0>  ?  R 1  <5  0>  ?  F4  <5  0) 

56'?  0=  1  ?  SO  '  5  0  1  ?  F' 1  5 O  '1  ?  S  E  T  H  <5  M  :•  ?  S" E T H G  v5  0 .■l  j  PS  I  C5  O.’1  ?  W  1 5  0  ?  YO  ( 5  OJ 

5700=  1 «  YETfiO  *  5  0>  ?  U  <  5  0>  ?  X  <50>  ?  CMPH  I  <5 0>  ?  SMPH I  <5 0> 

5710=  rx=co.  ?  i.:> 

5730=  Z=  CO.  ?  0.  > 

5730=  ZR=Z 

5740=  ZB=Z 

5750=  ZC=Z 

5760=  ZD=Z 

577  0=  MRXMM=MMMfi X - 1 

5780=  DO  3  MM= 1 ?  MRXMM 

5790=  M=MM-1 

5800=  R=3 

5810=  B=1 

5830=  C=1 

5830=  IFCMM.EQ. 1>  R=1 

584 0=  I F  <MM . EO . £>  B=£ 

585 0=  I F  <MM . EQ . 2>  C= 0 

5860=  10=IX^4-<-M> 

5870=  2=Z+R^YETR0  CMM> ♦CMPHI CMM> 

5380=  IFCMM.EQ. 1>  GO  TO  2 

589 0=  2R=2R+ I ©♦ <U < MM+ 1 ? ♦CMPH I CMM+ 1 > -B^U -MM- 1 > ♦CMPH I CMM- 1 > > 

5900=  1 ♦YO  CMM) 

59 1  0=  ZB=ZE+ 1  <U  <MM+ 1  ?  ♦SMPH  I  CMM+ 1 )  +U  CMM- 1  >  ♦SMPH I  CMM- 1  >  > 

5930=  1 ♦YO  CMM> 

593  0=  ZC=2C+ 1 £'♦  CX  CMM+ 1  >  ♦SMPH  I  CMM+ 1  >  -X  CMM- 1  >  ♦SMPH  I  CMM- 1  >  > 

594  0=  1 ♦YO  CMM> 

595  0=  ZD=ZD+ 1  ©♦  CX  CMM+ 1  >  ♦CMF'H  I  CMM+ 1  >  +C*X  CMM- 1  >  ♦CMF'H  I  CMM- 1  >  > 

5960=  1 ♦YO  CMM> 

5970=  GO  TO  1 

598  0=  £  Zfi=ZR+U  <2>  ♦CMF'H I  <£>  ♦YO  C 1  > 

599  0=  ZB=ZB+U C£> ♦SMPH I C2> ♦YO  C 1 > 

6  0  0  0=  ZC=ZC+X  <£>  ♦SMF'H  I  C3>  ♦YD  C 1  > 

6010=  ZD=ZD+X  C2>  ♦CMF'H  I  <£>  ♦YO  C 1  > 

6020=  1  CONTINUE 

6030=  3  CONTINUE 

6040=  RETURN 

6050=  END 
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SCATTERING 

CHQDGE 


IV  A  METALLIC  CIRCULAR  PICK 
—  VERSION  IS-'  17,-78> 


(TYPE  KA=0  TO  STOP  PROGRAM) 

CTVPE  KA=-1  FOP  P  DESPIPTION  OF  THE  PARAMETERS) 

(NORMAL I ZPT I  ON :  ESCAT=A^E I  NC^ENORM/  (£^R)  ♦EXP  |*K  ♦R-.  > 
CALL  ANGLES  IN  DEGREES) 


1. 

KA 

=  1 

c‘  • 

THETA  INCIDENT 

=45 

3. 

POLARIZATION 

=  0 

4. 

THETA  SCATTERED 

=45 

5. 

PHI  SCATTERED 

=  0 

WHICH  VARIABLE  IS  TO  EE  INCREMENTED?! 
TYPE  NUMBER  OF  CASES:  £5 
WHAT  IS  THE  INCREMENT?! 


KA 

CROSS 

SECTION 

E 

NORM 

SIGMA,-' 

<PI^A^£> 

THETA 

THETA 

PHI 

MAG 

PHASE 

1 

1. 00 

. 351E+00 

0. 

. 593E+  0  0 

-22. 65 

0. 

2.  00 

. 903E+00 

0. 

. 950E+00 

-110.86 

0. 

3.  00 

. 263E+00 

0. 

. 513E+00 

148.85 

0. 

4.  00 

. 353E+00 

0. 

. 594E+00 

85.  18 

0. 

5.  00 

. 334E+ 0 0 

0. 

. 578E+ 0 0 

-3.41 

0. 

6.  00 

. 447E+00 

0. 

. 669E+ 0 0 

-1 05.68 

0. 

7.  00 

. 149E+00 

0. 

. 387E+00 

148.80 

0. 

8.  00 

. 738E-01 

0. 

. 272E+00 

58.53 

0. 

9.  00 

. 659E-02 

0. 

. S12E-01 

43.29 

0. 

10.  00 

. 194E-01 

0. 

. 139E+00 

-45.88 

0. 

11. 00 

. 534E-01 

0. 

. 231E+00 

-95.93 

0. 

12.  00 

. 456E-01 

0. 

.21 3E+ 0 0 

167.71 

0. 

13.  00 

. 102E+00 

0. 

. 319E+00 

87.48 

0. 

14.  00 

. G27E- 0 1 

0. 

. 25  0E+ 0 0 

-11.15 

0. 

15.  00 

. 982E-01 

0. 

. 31 3E+00 

- 1 03 . £3 

0. 

16.  00 

. 637E-01 

0. 

. 252E+00 

159.30 

0. 

17.  00 

. 570E-01 

0. 

. 239E+ 0 0 

86.  09 

0. 

18.  00 

. 119E-01 

0. 

. 1 09E+00 

£1.45 

0. 

19.  00 

. 800E-02 

0. 

. S94E-01 

-127. 06 

0. 

£0.  00 

. 244E-01 

0. 

. 156E+00 

124.65 

0. 

21.00 

. 122E+00 

0. 

. 349E+  00 

88.27  • 

0. 

£2.  00 

. 172E+00 

0. 

. 414E+00 

71.90 

0. 

23.  00 

. 63 IE- 01 

0. 

. 251E+00 

1  0.  06 

0. 

£4.  00 

. 116E+00 

0. 

. 340E+00 

-120.81 

0. 

25.  00 

. 164E+00 

0. 

. 405E+00 

72.44 

0. 
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Y  AD-AOBO  367  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/6  17/9 

PREDICTED  MICROWAVE  ELECTRO-MAGNETIC  BACKSCATTERING  RETURNS  FOR— ETC<U> 
DEC  79  F  DX BARTOLOMEO 


UNCLASSIFIED  AFIT/GE/EE/79D-U  NL 


r~— 

s 

i 

■j 

1  00= 

PROGRAM  SCATTER  (  I NPIJT  •  OUTPUT-* 

'  i 

110= 

COMPLEX  A 1  •  A2 »  B 1 » B:£ » H  *  P 1  •  P2 .  F  •  U 1  *  U2 •  X*  Y»  ESCAT *  J 

; 

120= 

PERL  P I «  THETA*  A«  B*  D*  El  »  E£*  F 1  *  F£*  G1 » 62 *  Cl  *  C2 » KR» 

130= 

fS  1 » S£ »0»  T1 *  T2» RP6 *  V .  1.1 .  IEF'.  MMBSJ 0 *  MMBSJ 1 * MMBSJ2 

i 

1 4  0= 

I* I MENS I ON  ES CfiTM  •  1  0 0>  *  RNGLE  (  100>  *  S  I GMR  (100) 

i 

150= 

READ*  THETR 

} 

160= 

PI =3. 141593 

.| 

170= 

THETR=  (THETRX36Q.  > ♦ (£^PI> 

-! 

180= 

KR=0.  j 

1 

190= 

ft=l. 414214  j 

| 

200= 

£=-  <PI'"4> 

210= 

no  20  1  =  1  - £5 

j.1 

220= 

KR=KR+1 . 

230= 

£=£♦ < < <£^Kft> /P I> ♦♦ ( 1 . /£ . >  > 

; 

1 

240= 

E1=C0S < <PI/4. >  -  (THETR/2.>> 

, 

250= 

E£=COS  (  (P I /4 . >  +  <THETfl/2.>> 

1 

260= 

F1=RBS  <E1> 

| 

270= 

F£=RBS  (E£) 

280= 

G1=D^F1 

1 

, 

290= 

6£=D^F£ 

! 

300= 

CALL  FRESNEL  (G 1 » C 1  *  S 1  > 

310= 

CALL  FRESNEL CG2*  C£»  S2> 

320= 

H=CMPLX <0. » B>  ! 

' 

330= 

H=CEXP <H> 

340= 

P1*CMPLX<C1*S1> 

350= 

P£=CMPLX<C2»S£> 

360= 

R1=R^H^P1 

370= 

fi£=fl^H^P2 

380= 

RRG=2^KR» (SIN <THETR> > 

390= 

V=  (1  +  SIN  <THETR>  >  /  <2^SIN  (THETA.-)  > 

400= 

W= < 1  -  SIN  (THETA!)  > ✓ <£*S I N  (THETA>  > 

410= 

MMBS J2= <£/RRG> ♦MMBSJ 1 (RRG *  I ER>  -  MMBSJ 0 <RRG j I ER> 

'I 

420= 

X=  (  (A  1  ♦♦£>  ♦l.i.l  +  (R£^^£!)  ♦V>  ♦MMBSJ  1  (ARG*  I ER) 

1 

430= 

Y=  <  <R  1  ♦♦£>  ♦M  -  (fl£^£>  ♦V!)  ♦MMBS  J2 

i! 

44  0= 

J=  <0. » 1 . >  ! 

i\ 

i 

450= 

Y=Y  ♦  J 

460= 

Y1=RERL  <Y!) 

• 

470= 

Y£=R I MAG  <Y> 

480= 

X1=RERL  <X> 

490= 

X2=R I MAG  (X!> 

500= 

ESCAT1=X1-Y1 

510= 

ESCAT2-X2-Y2 

i 

520= 

ESCAT=CMPLX <ESCftT 1 * ESCAT2> 

530= 

ESC AT =ESCRT ♦ J 

1 

540= 

ESCAT =-ESC AT 

; 

550= 

ESCATR=REAL  (ESCAT> 

i 

56  0= 

ESCAT  I  =AIMAG  (ESCAT!) 

i 

570= 

ESCATM  ( I >  =  (  (E SCATP^^£ .  >  +  (ESCAT I ♦♦£.>>♦♦ <1 ./£. > 

580= 

SIGMA <I>=  (ESCATM ( I >♦♦£.!> 

590= 

ANGLE=ATAN  (ESCAT  I -'ESCRTR} 

600= 

ANGLE (I> =( (ANGLE) ^(£^PI>> *360.  : 

610= 

IF < (ESCATP.LT. 0. > . AND. (ESCAT 1 . GT. 0. >  >  ANGLE <I> =ANGLE CD  +  V 

80. 

• 

620= 

IF((ESCATR.LT. 0.> .AND.  (ESCATI.LT. 0.>>  ANGLE  CD  =  ANGLE  <D  -  1 

CO 

o 

• 

1 

630=  20 

CONTINUE 

640= 

KA=0. 

650= 

PRINT30 
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6*0=  30  FORMAT  <11 X >  // » 8X »  “KA“»7X»  14X»  “ESCRT"*/* 

67  0=  $37X » " MAG “ . 9X « “ PHASE " . ✓> 

680=  DO  4  0  1  =  1 j  £5 

£ Q  |"|=  ^  ^ 

7  0  6=  PP I NTS 0 » K  A «  S  I GMA  < I > » ESC  ATM < I > , ANGLE  < I > 

710=  SO  FORMAT  <"  " , 3X. F7. 2- 5X« El  0. 3, 7X» El  0. 3f 3X» F7. 2> 

780=  40  CONTINUE 

730=  STOP 

740=  END 

750=  SUBROUTINE  FRESNEL <AA. C» S> 

760=  D I  MENS  I  ON  X  <2  0  0  0>  .  Y  <£  0  0  0>  •  ZX  <2  0  0  0>  *  ZY  <2  0  0  0> 

770=  INTEGER  NDIM 

780=  IF  < CAR. 6T. -0. 01) . AND. <RA. LT. 0. 01>> GD  TD  6 

790=  GO  TO  11 

800=  6  C=0.  0 

910=  S=0.  0 

S£0=  GO  TO  4 

830=  11  H=. 01 

84  0=  BB=ABS  <AA> 

850=  NDIM=BB/H  +  .1 

860=  K=NDIM+1 

870=  DO  10  1=1 »K 

880=  A=I^H  -  H 

890=  X  <  I  >  =COS  <  1 . 57  08  ©♦  <  A^2 .  >  > 

9 0  0=  Y  < I >  =S I N  < 1 . 57  08  0*  < ^♦♦2 . >  > 

910=  10  CONTINUE 

920=  SUM2-0. 

930=  SUM4=  0 . 

940=  1  HH=.5^H 

95  0=C 

96 0=C  INTEGRATION  LOOP 

970=  DO  2  1=1 » NDIM 

980=  SUM 1 =SUM2 

990=  SUM3=SIJM4 

1000=  SUM2=SUM2  +  HH»  <X < I > +X  < I  + 1 >  > 

1010=  SUM4=SUM4  +  HH^<Y <I>+Y <I+1>> 

1020=  ZX  <I>  =SUM1 

1030=  2  Z Y  < I >  =SUM3 

1 04  0=  3  ZX  <ND I M> =SUM2 

1050=  ZY  <NDIM>  =SUM4 

1060=  7  IF  <  AA. GE. 0. 1> GO  TO  8 

107  0=  C=  <- 1 .  >  ♦ZX  <ND  I  M> 

1 08 0=  S=  <- 1 . > ♦ZY  <NB I M> 

1090=  GO  TO  9 

1100=  8  C=ZX  <NDIM> 

1110=  S=ZY  <NDIM> 

1120=  9  CONTINUE 

1130=  4  RETURN 

1140=  END 
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j 

! 


Kfl  S I GMfi 


ESCfiT 

MRG  PHRSE 


1.00 

2.  00 

3.  00 

4.  00 

5.  00 

6.  00 

7.  00 

8.  00 
9.  00 

10.  00 
11.00 
12.  00 

13.  00 

14.  00 

15.  00 

16.  00 

17.  00 

18.  00 

19.  00 

20.  00 
21. 00 
22.  00 
23.  00 
24.00 
25.  00 


. 409E+00 
. 104E+01 
. 370E+00 
. 700E+00 
.  297E+  0  0 
. 265E+00 
. 848E-01 
. 714E-01 
. 226E-01 
.240E-01 
. 444E-01 
. 340E-01 
.980E-01 
. 652E-01 
. 110E+00 
. 563E-01 
. 676E-01 
. 253E-01 
. 275E-01 
. 169E-01 
. 138E-01 
. 294 E- 01 
. 199E-01 
. 513E-G1 
.330E-01 


.  640E+00 
.  1 02E+01 
.  608E+00 
.  837E+00 
.  545E+00 
.  514E+00 
.  291E+00 
.  267E+00 
.  150E+00 
.  155E+00 
.  21 1E+00 
. 1S4E+00 
. 313E+00 
. 255E+00 
. 332E+00 
. 237E+00 
. 26  0E+ 0 0 
. 159E+00 
. 166E+0O 
. 130E+00 
. 117E+00 
. 171E+00 
. 141E+00 
. 226E+00 
. 182E+00 


-.09 
-92.45 
2.59 
105.75 
-165.44 
-67. 00 
14. 66 
100.69 
147.33 
-120. 14 
-80.24 
12.  08 
90.71 
-178.69 
-85.81 
1.94 
94.77 
168.41 
-99. 14 
-51.91 
45.27 
97.22 
-174.10 
-91.91 
-5.  07 


! 

i 

j 

1 

j 


i 


| 


) 


i 


i 


i 
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I  0 0=  PROGRAM  P I  SK£5  ( I NPUT ,  OUTPUT.  TAPE5=  I NPUT »  T APE6*0UTPUT> 

I I  6=c  FAR  FIELD  SCATTERING  BY  A  CIRCULAR  METALLIC  DISK 

1£0=  DIMENSION  LABEL  C5>  *  VAR  <5> 

1 3 0=  COMPLEX  F4 , PS I , M . VO , YETAO *  U «  X ,  Z «  ZA ,  ZB ,  ZC  » ZD 

140=  1,1  X ,  EPART ,  EF'APP .  EPEPT ,  EPEPP ,  ESNT ,  ESNP 

1 5 0=  COMMON  EI6  <5 CO  ,  D  C5 0 ,  5 CO  ,  DNEG  1  5 0>  ,  P 1  <5 CO  ,  F4  C5 0> 

160=  l»SOC50>  ,PC50>  »SETAC50>  ,  SETAO  <5  CO  >PSI  C50>  »W<50>  »YO<50> 


170= 

1 ,  YETAO  <5  0  >  ,  U  5  CO  ,  5 

C5  CO  ,  CMPH I  C5  0>  ,  SMF'H  I  C5  CO 

180= 

DATA  LAPEL  Cl > -a  OH 

KA 

/ 

190= 

DATA  LABEL  C£>  s 1  OH 

THETA  I 

/ 

200= 

DATA  LABEL  C3>  a  OH 

POL 

/ 

210= 

DATA  LABEL C4)/l OH 

THETA  S 

/ 

220= 

DATA  LABEL C5>/1 OH 

PHI  S 

/ 

230= 

XX* (0. , 1.) 

24  0= 

WRITE  C6,  U> 

25  0= 

11 

FORMAT  C 1 X , s// , IX,  1 

'SCATTERING  BY  A  METALLIC  CIRCULAR 

260= 

I/'jSX,  "  CHODGE  —  VERSION  12 

yl7y7S) "> 

270= 

WR I TE  C6 *  26) 

280= 

26 

FORMAT  Cl Xj-w,  IX, ' 

"  CTYPE  KA 

=0  TO  STOP  PROGRAM)",/, 

£90=  1 1 X ,  "  CTYF'E  KR=- 1  FOP  fi  DESK  I PT I  ON  DF  THE-", 

3  0  0=  1 "  PARAMETERS'1  " ,  / ,  1 X ,  “  (NORMAL  I ZAT I  ON :  ESCAT=A»E  I  NC^ENOPf ' 

31  n=  1  ‘V  C£^R>  ♦EXP  C— -.•♦K ♦R)  >•*,/,  IX,  “  CALL  ANGLES  IN  DEGREES.-'  "> 

320=  4  NOC= 1 
330=  I NDEX= 1 

340=  WRITE  <6, 27> 

35  0=  27  FORMAT  < 1 X ,  /v  ✓ , 1 X , " 1 .  KA" , 14X, “=  “ > 

36 0=  READ <5 , ♦>  V AR  < 1 > 

370=  IF  <VAR  Cl>  . EQ. -1>  GO  TO  41 

38  0=  I F  < VAR  <1 > . LE . 0 . > 60  TO  5 

390=  WRITE C6,£8> 

400=  28  FORMAT (IX, "2.  THETA  INCIDENT  =  "> 

4 1  0=  READ  C5 ,  ♦>  VAR  <2> 

420=  WP I TE  •- 6 , 29) 

430=  29  FORMAT <1X, "3.  POLARIZATION  =  "> 

44  0=  READ <5 , ♦>  VAR <3> 

450=  WRITE  <6, 3  0> 

460=  30  FORMAT  «C1X,  ”4.  THETA  SCATTERED  =  "> 

470=  READ  <5,  ♦>  VAR  C4> 

480=  WRITE <6, 31 > 

490=  31  FORMAT  ax,  "5.  PHI  SCATTERED  =  "> 

5  0  0=  READ  <5, ♦>  VAR  C5> 

510=  WRITE C6,3> 

520=  3  FORMAT  (IX,//, IX, "WHICH  VARIABLE  IS  TO  BE  INCREMENTED? " > 
530=  READ  <5, ♦>  NVAR 

540=  I F  C  CNVAR .  LE .  CO  .  OR .  (NVAR .  GT .  50  >  GO  TO  23 

550=  WRITE <6, 21 > 

560=  21  FORMAT C1X, "TYPE  NUMBER  OF  CASES: "> 

570=  READ  C5,’  ♦>  NOC 

580=  WRITE  <6, 32> 

590=  32  FORMAT <1X, "WHAT  IS  THE  INCREMENT?") 

600=  READ  <5» ♦)  VINCRE 

610=  23  CONTINUE 

62  0=  WR I TE  (6 , 24)  LABEL  CNVAR) 

630=  24  FORMAT  <  1 X ,  /y ,  A 1  0 , 4X ,  "CROSS  SECTION"  ,£2X,  "E  NORM" ,  ✓, 

640=  11 3X ,  "SI GMA/-  CP  I  ♦A^£)  " ,  1 1 X ,  " THETA " .  1 5X ,  "PHI",-'',  1 3X , 

650*  1 "THETA" , 7X, "PHI",8X, ”MAG",5X, " PHASE ">6X, "MAG",5X, 


**  0® 
*70=  * 
*80= 
*90= 
700= 

71  0= 
720= 
730= 
740= 
750= 
7*0= 
770= 
780= 
790=  34 
800= 
810= 
820=  ’ 
830= 

84  0= 
850= 
8*0= 
870= 
880= 
890= 
900= 
910= 
920=  7 
930= 

94  0= 
950= 
9*0= 
970=  10 
980= 
990= 
1000= 
1010= 
1020=  13 
1030= 
1040= 

1  050= 
10*0= 
1070= 
1080= 
1090= 
1100= 
1110= 
1120=  2 
1130= 
1140=  8 
1150= 
11*0= 
1170= 
1180= 
1190= 
1200= 


1 "PHASE",/) 
c=var  <:i> 

NMMAX=45 

NNMRX=45 

IRPMRX=45 

THEO=VRP  <£)  ♦S.  14159 ''ISO 
ETRD=COS (THED) 

THE=VAR  <4)  ♦S.  14159/180 
ETR=CDS (THE) 

PH I =VAP <5> ♦3.141 59/ 1 8  0 
RLF=VAR C3> ♦S. 14159/180 
CALF=COS <ALF) 

SRLF=SIM  <RLF> 

DO  10  MM= 1 ,  MMMRX 

M=MM-1 

IQ=0 

CRLL  SMNO  <M , MMMRX) 

CRLL  DEE I GM < C  ,  M , MMMRX) 

CRLL  DBCOFN <C • M , MMMRX . MMMRX , I RPMRX> 

CRLL  DNEGN CC »  M » MMMRX) 

CRLL  OBRRD <C , M» IQ* MMMRX  *  MMMRX , I RRMRX)  - 
IF  <10.  EC!.  l)GO  TD  10 
CRLL  FPSI CM. MMMRX) 

CRLL  PDLYM CETRO » M» I RRMRX) 

CRLL  OBPNG <MMMPX , I RRMRX) 

BO  7  1=1, MMMRX 
SETRO  < I ) =SETfl  < I ) 

CRLL  FI,. I  CM,  MMMRX) 

CALL  PDLYM (ETA, M» I RRMRX) 

CRLL  OEPNG (MMMRX » I RRMRX) 

CRLL  FY(M, MMMRX) 

CONTINUE 

DO  13  MM=1, MMMRX 
CRLL  FXUvMM, MMMRX) 

IF (MMMRX. LT. MM) GO  TD  13 
CRLL  CSPHX<MM»PHI> 

CONTINUE 

CRLL  FZ  (MMMRX , Z , ZR, ZE, ZC, ZD) 

EPRRT=ETR4  <-£+Z+CHPHl <£) +ZR) 

EPRRP— £^Z^SMPHI  <£)  +ZE 
EPERT=ETfi^  <£^Z*SMPH I <2) ~ZC) 
EPERP=-2*Z*CMPHI <£) +ZD 
IF(ETRD.EQ.  0.)60  TD  2 

ESMT=2^IX^ CCRLF^EFRRT/ETRD+SRLF^EPERT) /C 
ESNP=£»IX*  ’.CRLF^EPRRR/ETRD+SRLF^EPERP)  /C 
GD  TD  8 

ESMT =£♦ I X^SRLF ♦EPEPT /C 
ESMP=£^IX^SRLF*EPERP/C 
EMAGT=CRBS (ESMT) 

EMRGP=CABS (ESNP) 

S I GTHE=EMRGT ♦EMRGT 
S I GPH I =EMAGP^EMR6P 
E1=RERL(ESNT> 

E£=RIMRG (ESMT) • 

IF  (El .  EG!.  0. )  GD  TD  16 
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1210=  ARG=E£^E1 

1  ££ 0=  EPHAT  =  1 8  0--  3 .141 59^ATAN  CAPG> 

1 £3 0=  I F  C CE 1 . LT . 0 . > . PHD . CE£ . GT . 0 . > >  EPHAT =EPHAT + 1 8  0 

1  £4  0=  I F  C  CE  1 .  LT .  0.  >  .  AND .  •  ££.  LT .  0 .  >  >  EPHAT=EPHAT- 1 8 0 

1 £5  0=  GO  TO  17 

1£60=  16  EPHAT =90 

1  £7 0=  IFcE£.LT.  0. > EPHAT=-90 

1 £3  0=  17  E1=REALCESNP> 

1 £9  0-  E£=A I MAG  CESNP> 

1300=  IF  <E1 . EQ. 0. >GD  TO  18 

1310=  AR6=E£--E1 

13£0=  EPHAP= ISO'S.  14159^ATAN  CARG> 

1330=  IFCCEl.LT.  0.>  .AND.  CE£ .  GT .  0 .  >  >  EF’HRP=EPHRP+ 180 

1340=  IFCCE1.LT. 0.> .AND. CE£.LT. 0. > > EPHRP=EPHRP-180 

1350=  GO  TO  19 

1360=  18  EPHAP=90 

1370=  IFCE8.LT. 0. > EPHRP=-90 

1 38  0=  1 9  IF  CNVAR . EQ . 0> NVRP= 1 

1 39  0=  I F  CVAR C 1 > . GT . 1 5>  GO  TO  1000 

1400=  GO  TO  £000 

1410=  1000  CALL  SCATTER CVAR- EM AGT* EPHAT* SIGTHE> 

14£0=  £000  CONTINUE 

1 43 0=  WR I TE  C6 - £5>  VRR CNVAR> *  S I GTHE - S I GPH I > EMRGT * EPHRT « 

1440=  1EMAGP-EFHAP 

1450=  £5  FORMAT C1X* F7. £- £ C1X* El  0. 3>  -£C1X* El 0.3* 1X*F7.£>> 

1460=  37  IFCINDEX.EQ.NOOGO  TO  4 

1470=  INDEX-INDEX+1 

1 48  0=  VAR  CNVAR> = V AR  CNVAR>  + V I NCR E 

1490=  60  TO  6 

1500=  5  CALL  EXIT 

1510=  41  WRITE  C6*  40> 

1 5£ 0=  40  FORMAT  ClX*/v*  IX *  "THE  DISK  DF  RADIUS  A  LIES  IN  THE  X-Y  PLANE 

1530=  1 IX* "CENTERED  AT  THE  ORIGIN.  THE  CONVENTIONAL  CR*  THETA*”*, 

IX* 

1540=  1  “  PH  I  >  COORDINATE  SYSTEM  IS  USED  IN  THE  FAR  FIELD.  "j/'jIX* 

1550=  1 "THE  PLANE  OF  INCIDENCE  DF  THE  PLANE  WAVE  IS  THE" *  •'j IX* 

1560=  1"X-Z  CPHI=0>  PLANE.  THE  POLARIZATION  ANGLE  CPOL> " * .'*  IX- 

1570=  1  "OF  EINC  IS  MEASURED  FPDM  THE  PLANE  OF  INCIDENCE  "*.'*1X* 

1580=  1 "  I N  THE  PH  I -DIRECT  I  ON.  I.E.*  PDL=0  IS  THE  PARALLEL"  *  •'*  IX* 

1590=  1 " CTHETA>  CASE  AND  POL=90  IS  THE  PERPENDICULAR  CPHI) "*✓* IX* 

1600=  1"CASE. " *£X* "THE  RESULT  IS  A  SOLUTION  OF  THE  RIGOROUS" *✓> 
1610=  1 1 X  * " E I GENFUNCT I ON  SCATTER  I N6  PROBLEM . " > 

16£0=  GO  TO  4 

1630=  4£  CONTINUE 

1640=  GO  TO  4 

1650=  END 

1660=C 

1670=C  OBLATE  SPHEROIDAL  ANGULAR  FUNCTION* S* OF  ARGUMENT  0* 

1680=C  ORDERM-N!  WITH  N-M  EVEN* UP  TO  ORDER  N=M+£^NNMAX-£. 

1690=C  CEQUAL  TO  PMNC0>> 

1700=C 

1710=  SUBROUT  I NE  SMNO-CM  *  NNMAX> 

17£0=  COMPLEX  F4 

1 73  0=  COMMON  E I G  C5  0>  *  D  C5  0  *  5  0>  *  DNEG  C5  0)  *  R 1  C5  0>  *  F4  C5  0> 

1740=  1  *  SO  C50> 

1750=  S0C1>  =  1 
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1760= 

1770= 

17:30= 

1790= 

1 8  0  0= 
1810= 
1880= 
1880= 

1 84  0=< 

1 85  0=1 
1880=1 
1870= 
1880= 
1890= 
1900= 
191  0= 
198  0= 
1930= 
194  0= 
1950= 
1960= 
1970=1 
1980=1 
1990=< 
£000=1 
8010= 
8080= 
8030= 
8040= 
8050= 
8060= 
8070= 
8080= 
£090= 
8100= 
8110= 
8180= 
2130= 
2140= 
8150= 
2160= 
8170= 
2180= 
2190= 
8800= 
2210= 
2880= 
2830= 
2840= 
2850= 
8860= 
8270= 
2880= 
8890= 
2300= 


1 


IF  r'M.EO.  0.>GO  TO  1 
DO  £  MM= 1  •  M 

£  SO  1  '>  =  • £^MM- 1  >  ♦SO  <  1  > 

1  DO  3  NN=  1 ,  NNMAX 
N=£^<NN-1>+M 

9  30  <NN+ 1  =-  CN+M+l  >  »SO  CNN}  ✓  <N-M+2> 

RETURN 

END 

SIN  AND  COSINE  FUNCTIONS  OF  PHI 

SUBROUT I NE  CSPH I <MM •  PH  I  > 

COMPL  EX  F4  *  PS  I  *  I-' •  YO  >  YETfiD  •  U  j  X 

COMMON  EIG  f«!n>  *  D  <’50»  50>  •  DNEG  <50>  •  R1  <50>  » F4  <50> 

1 ,  SO  <  5  0>  *  P  <5  CO  ,  SETA  <5  0>  » SETfiO  <5  0>  » PS  I  <50>  » U  <5  0>  » YO  <5  0> 

1  *  VETAO  '"5  0>  ?  U  <5 0>  *  X  <5 0>  » CMPH I  <5 CO  •  SMF'HI  <5  0> 

M=MM— 1 

CMPH I <MM> =COS CM^PHI) 

SMF'HI  <MM>  =S I N  -  M^PHI) 

RETURN 

END 

OBLATE  SPHEROIDAL  EIGENVALUES  DF  ARGUMENT  C*  ORDER  M»N 
WITH  N-M  EVEN  UP  TO  ORDER  N=M+2^NNMAX-2 


4 


2 


3 


13 


SUBPOUT I NE  QBE  I GN  <C *  M < NNMAX> 

COMMON  E I G  <5  0> 

DIMENSION  IP  <50>  *  ALPHA  <50>  >  BETA  <50>  » P  <50> 

CONTINUE 

M8=8^M 

C2=C^C 

ACC= 1 . 0E-05 

NN£=NNMAX+£ 

N I =NN2+ 1 
P<1>=1 


IPC1>=1 

DO  £  IGO=l»NN£ 

I  *•/=£♦  I  GO- 1 

IW=M£+£^IGQ 

IX=M£+4MG0-1 

ALPHA  <  I Q CO  =  <C£+  <M£^  <£♦  I V- 1  >  +£♦  I V*  <  I V- 1  >  - 1  >  >  /'  <  I X^  <  I X-4>  > 
t- CM+IV-1> ♦ (IV+M) 

BETA  <  I QG+ 1  >  =C2^  I X^SGRT  <  I ’•.•’♦  <  I V+ 1  >  ♦  I ••)♦  <  I W- 1  >  /  C I X^  I X-4 .  0>  > 
BETA  <NN8+ 1 >  =  0 . 

EO=ABS < ALPHA <1>  >  +ABS <EETA <£>  > 

DO  3  IGG=£»NN£ 

AD=AES CBETA < I GG>  > +AES C ALPHA < I GG> >  +ABS <BETA < I 00+ 1>> 

BETA  *:'IGG>  =BETA  <IGG>  ♦BETA  CIOGO 
IF  CRO.GT. BO>  BO=AO 


CONTINUE 

AO=-BO 

BDI=BO 

CONTINUE 

BO=BOI 

DO  20  I GQ= 1 » NNMAX 
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£310= 

N=£*IOO-£+M 

£3£  0= 

A=AD 

£33  fi= 

B=BO 

£34  0= 

IEPR=-1 

£35  0= 

£1 

IIS=0 

236  0= 

CO=  ‘  A+B>  ■'■  £ 

£37  0= 

I F  <CO>  5  0  *  ££  *  5  0 

£380= 

50 

ERR=  CB-A>  .-'AES  <CO> 

£390= 

IERR=IERP+1 

£400= 

IF  < I ERR-6 0>  40*41*41 

£41  0= 

41 

WRITE  (6*  4£>  N 

£420= 

4£ 

FORMAT C£X* "ITERATIONS  EXCEEDED  FOR  EIGENVALUE  "*I3> 

£430= 

GO  TO  700 

£44  0= 

40 

I F  CERR-ACO  £4  *  £4 ,  ££ 

£450= 

££ 

P  <2> =ALPHA < 1 > -CO 

£460= 

DO  5  1=3* NI 

£470= 

P  <I>  =  (ALPHA  I  —  1  >  -CO-BETA  <I-1>  ♦  <P  CI-£>  /P  <I-1>  >  >  <I-1)  ! 

£48  0= 

PMAG=ABS<P<I>>  i 

£490= 

IFCPMAG.GT. 1. 0E+33>  GO  TO  7  ! 

£500=  5  CONTINUE 

£510=  1£  CONTINUE 

£5£0=  no  6  I =£  >  N I 

£530=  IF  <Pa>>  14*8*9 

£540=  $  IF <P<I-1>> 9*9*14 


£550= 

14 

IP<I>=-1 

£560= 

GO  TO  10 

£570= 

9 

ipa>=i 

£580= 

10 

IF(IPa>-IP<I-l>>6*  11*6 

£590= 

11 

IIS=IIS+1 

£600= 

6 

CONTINUE 

£610= 

IFaiS-IGCO  16*  15*  15 

2620= 

1 5 

a=co 

£63  0= 

GO  TO  21 

£64  0= 

16 

E=CO 

£650= 

GO  TO  £1 

2660= 

£4 

BO=CO 

2670= 

700 

EIG  < I GQ>  =-CO 

£680= 

£0 

CONTINUE 

2690= 

RETURN 

£700= 

7 

NNMAX* 1—4 

2710= 

N I =NNMAX+3 

2720= 

GO  TD  4 

£730= 

END 

£740=C 

£750=C  OBLATE  SPHEROIDAL  EIGENFUNCTION  EXPANSION  COEFFICIENTS 
£760=C  OF  ARGUMENT  C*  ORDER  M  *  N  ?  R  ?  WITH  N-M  EVEN*  UP  TO  ORDER 

£77  0=C  N = M + £ ♦ N N M A X - £ *  P=£^ I RPMAX+£ 

£780=C 

279  0=  SUBROUT I NE  OBCOFN  <C  *  M  *  MMMAX » NNMAX *  I RRMAX> 

£800=  COMMON  EIG  <50>  *  D  <T50*  50) 

28 10=  D I  MENS  I  ON  DP  <5  0> 

2820=  C£=C^C 

2830=  MM=M+1 

2840=  DO  1  NN=1* NNMAX 

2850=  N=M+£«-NN-£ 


2360= 

4 

DPCIPPMRX+3>=0 

2870= 

HP  C I RPMRX+2)  =  1 . 0E-3 0 

2880= 

d  cnn>  n=o 

£0Qfl= 

D  CNN >  2)  =  1 

2900= 

jj=  cn-m)  /£+i 

2910= 

I'D  107  LL=1  >  IPPMRX 

2920= 

L=LL-1 

2930= 

IF  CLL. GE. JJ>  L= I PRMRX+ J J-LL 

294  0= 

IR=2H. 

2950= 

IRM=n+IR 

2980= 

RR=  < M+ 1  RM+2)  ♦  *  ri+ 1 RM+ 1  >  ♦CE/  C  <£♦  I  RM+3)  ♦  <£♦  I  RM+5)  ) 

2970= 

ER= <2*IRM« C I RM+ 1 > -£^M^M-1) ♦Cc'/  C  C2MRM-1)  ♦ 

2980= 

1  <£♦  I  RM+3)  )  - 1 RM+  C  I PM+ 1  ) 

2990= 

CR=IR+  C I R—  1 )  +C£/  C  <£♦  I RM— 3)  ♦  C2+IRM-1) ) 

3000= 

IF  CLL— J  J)  1  05  >  1  08  >  1  08 

3010= 

105 

D  CNN « L+3> =- CCR+D  CNN > L+ 1 )  +  CEP+E I G CNN) ) +D CNN >  L+2) 

3020= 

DMRG=RBS CD  CNN >  L+3) ) 

3030= 

IFCDMRG.GT. 1. 0E+30)  GO  TO  3 

304  0= 

GO  TO  107 

3050= 

108 

DP CL+1) =- CRR+DP CL +3)  +  CER+EIG  CNN) > ♦DP  CL+2)  ) /CR 

3080= 

DMRG=RES  CDP  CL+l) > 

3070= 

IF  CDMRG. GT. 1 . 0E+30)  GO  TO  3 

3080= 

107 

CONTINUE 

3090= 

DL=RBSCDCNNf JJ+1)> 

3100= 

DL=RL0G1 0  CDL) 

3110= 

DLP=RES  CDP  C.JJ+1)  ) 

3120= 

DLP=RL0G1 0 CDLP) 

3130= 

DL=RES CDL) 

3140= 

DLP=RES  CDLP) 

3150= 

DL=DL+DLF 

3180= 

IF CDL. GT. 30.)  GO  TD  5 

3170= 

CON=D CNN > J J+ 1 ) /DP C J J+ 1 ) 

3180= 

RCON=RES CCON) 

3190= 

IFCRCON.GE. 1. OE+32)  GO  TD  2 

3200= 

DO  118  J=JJ> IRRMRX 

3210= 

118 

D  CNN  j  J+2) =CON+DP CJ+2) 

3280= 

F=1 

3230= 

IF  CM) 198? 198 > 199 

3240= 

199 

DO  110  1  =  1  >M 

3250= 

110 

F=F+CM+I) 

3280= 

198 

SIJM=  0 

3270= 

MMX= I RRMRX+ 1 

3280= 

DO  113  1  =  1 j  MMX 

3290= 

IR=2+I 

3300= 

SUM=SUM+F+D CNN  > 1  +  1) 

3310= 

IFCI-JJ)  11 3 >19? >113 

3320= 

197 

FNM=F 

3330= 

113 

F=  C-F^ C I P+2+M- 1 ) ) / 1 R 

3340= 

RLF=FNM/SUM 

3350= 

DO  114  I  =  1 >  MMX 

3380= 

D  CNN  > I) =RLF+D  CNN> 1  +  1) 

3370= 

114 

CONTINUE 

3380= 

1 

CONTINUE 

3390= 

RETURN 

3400= 

3 

IRRMRX=LL-1 

3410=  GO  TO  4 

342  0=  £  IF'PMAX=IPRMRX-1 

3430=  GO  T04 

344  0=  5  NNM AX=NN- 1 

3450=  RETURN 

346  0=  END 

347 0=C 

348 0=C  NEGATIVE  D  COEFFICIENT  SUBROUTINE 

34Rft=C 

35  0  0=  SUE POUT  I ME  PHEGN <C  *  M » NNMAX1 

35 1  0=  COMMON  E I G  <5  01  -  P  <5  0-501- PNEG  <5 0> 

35£0=  PO  4  NN=1 » NNMAX 

3530=  C£=C^C 

3540=  IF  CM. GE.  11  GO  TO  £ 

3550=  PO  5  1 1  =  1 »  NNMfiX 

358  0=  5  PNEG  C 1 1 1 =P < 1 1 - 1 > 

3570=  GO  TD  3 

3580=  £  B 1=1.0 

3530=  B£=0. 0 

3600=  EI=EIGCNN1 

3610=  PO  1  IRR=1 > M 

36£  0=  IR=£#IRR-S#M-S 

3630=  AR= (£^M+IR+£> ♦ C2*t!+IR+11 ♦C£/ < <£*M+£^IR+3> 

364  0=  1 ♦ C£^M+£M R+51  > 

365 0=  BR= <M+ I R> ♦ • M+ IR+11 -El- <£♦ C M+ I P> ♦ <M+ I P+ 1 > -£^M*M- 1 1 

3660=  1  ♦C£/  <  <£^M+£-»  I R- 1  >  ♦  <£^M+£^  I  R+3>  > 

367  0=  CR= C I R> ♦ C I R- 1 1 ♦C£/ C  <£>M+£^I R-3> ♦ C£^M+£^ I R- 1 1 > 

3680=  B3=B£ 

3690=  E£=E 1 

3700=  1  Bl=-:EP^B£-CP>B3>/fiR 

3710=  R=P<NN»  11 -'Bl 

37£0=  PNEG  CNN1  =ft 

373 0=  PPM=ABS  '-PNEG  CNN1  > 

3740=  IFCPPM.LT. 1. 0E-35>  GD  TO  6 

3750=  4  CONTINUE 

3760=  3  RETURN 

3770=  6  NNMAX=NN-1 

3780=  RETURN 

3790=  ENP 

3800=C 

381 0=C  OBLRTE  SPHEROIPRL  RRPIRL  FUNCTION  R  <4>  OF  RPGUMENT 

38£0=C  ORPER  M-N?  WITH  N-M  EVEN?  UP  TO  ORPEP  N = M + £ ♦ N N M R X - 

383 0=C  ALSO  NORMALIZATION  FUNCTION.-  N. 

384  0=C 

385  0=  SUBROUT I NE  DBRAP <C - M -  10- MMMRX - NNMRX > I PPM AX> 

386  0=  COMMON  E I G  C5  0  •  >  P  <5 0  -  5  0>  -  PNEG  <5  01  -  R 1  <5  0>  -  F4  <  5  01 

387  0=  COMPLEX  I X - R4 - F4 

3880=  IX=C0. 0-1.01 

3890=  EERC=1 

3900=  EAC=1 

3910=  EAC£=1 

3920=  GRO=l 

3930=  IF  CM. EQ. 01  GO  TO  £0 

3940=  MRXM=M+1 

3950=  PO  19  MM=£- MAXM 


3960= 

IM=MM— 1 

39?  0= 

GR 0=  ■  £♦  I N- 1  )  ♦  •:  £♦  I  M>  ♦GP'D 

3980= 

EERC=  <2»I U- 1 > ♦EE6C 

399  0= 

EfiC2=  <£♦  I  r1- 1  >  ♦  <£+ 1  M)  ♦ERC2 

4000= 

19 

EhC= I N^EhC 

4010= 

I F  CEERC  .  GE .  1 .  0E+ 1  ?>  GO  TO  4 

4020= 

IF cERCS. GE. 1 . OE+30)  GO  TO  4 

4  030= 

20 

DO  17  NN=1, NNMRX 

4  040= 

N=24 CNN-1) +M 

4  05  0= 

SUM=  0 

4  06  0= 

RPsRPQ 

4070= 

ENOPM=0. 

4080= 

I'D  18  NR=1, IRRMRX 

4  090= 

IR=2^CNR-1) 

4100= 

SUMP=GR^DCNN»NR) 

4 1 1  0= 

SUM=$UM+SUMP 

4120= 

DMRG=RES  CD  CNN, NR) > 

4130= 

IFCDNRG.LT.  1.  0E-30)  GO  TO  1 

414  0= 

EN0PMP=2-*6R^  CD  CNN ,  NR)  )  ♦♦2--''  <!£’♦  I R  »-2+M+ 1 ) 

4150= 

CL 

ENDPM=ENDRM+ENORMP 

4160= 

G  R  M  R  G = R  B  S  C  G  R  ) 

4170= 

IF CGRMRG. GT. 1 . 0E+30)  GO  TD  5 

4180= 

GR=  I  R+2^M+ 1 .  )  s  C I R+ 1 . )  ♦  C I R+2+N+2 .  )  /  C I R+2 .  )  ♦GR 

4190= 

18 

CONTINUE 

4200= 

R 1  CNN>  =  C- 1  >  ♦♦  CNN- 1  >  ♦^•♦♦M^ERC^C^M+D  CNN ,  1  >  /  <  C2^M+ 1  )  ♦SUM) 

4210= 

P2=  C- 1  >  ♦♦  CNN- 1  >  ♦  C2^M- 1  >  ♦ERC^C^  CM- 1  >  ✓  Cc’^ERCE)  4-3 .141 59 

4220= 

!♦  CEERC^E'SIJM)  ^♦♦M/'DNEG  CNN> 

423  0= 

RR2=RBS CR2> 

424  0= 

IFCRR2.GE.  1. 0E+30)  GO  TD  4 

4250= 

EERC= CN+M+ 1 ) ♦EERCv  CN-M+2) 

4260= 

I F  CEERC . GE . 1 . 0E+ 1 7>  NNMRX=NN 

4270= 

R4=R 1  CNN) -IX4R2 

4280= 

R5=RBS CENORM) 

4290= 

R1=RL0G1 0  CR5> 

4300= 

R4=PBS CR2> 

4310= 

R2-RL0G1 0 CR4) 

4320= 

R3=R1+R2 

4330= 

IFCR3.GT.30.)  GO  TO  3 

434  0= 

R5=RBS  CR1 CNN)) 

4350= 

fl 1 =RLOG 1 0 CR5) 

4360= 

R3=RBS Cfi 1 ) +RES CR2) 

4370= 

I F  CR3 . GT . 3 0 . )  GO  TO  3 

4380= 

F4  CNN)  =  1  -c  CENORM4R4) 

439  0= 

17 

CONTINUE 

4400= 

RETURN 

4410= 

1 

ENOPMP=0 

4420= 

GO  TO  2 

443  0= 

3 

F4  CNN)  =  CO. , 0.) 

444  0= 

GO  TO  17 

4450= 

5 

GP=0 

4460= 

GO  TO  18 

4470= 

4 

MMMRX=M 

443  0= 

IQ=1 

4490= 

RETURN 

4500= 

END 
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451 0=C 

452  0=C 

453  0=C 

454  0= 
4550= 
456  0= 
4570= 
45?  0= 
4590= 
4600= 
4610= 
4620= 
463  0= 
4640=  1 

465  0= 

466  0= 
467 0=C 
468  0=C 
469 0=C 
4700=0 
4710= 
4720= 
473  0= 
4740= 
4750= 
4760= 
4770= 
4780= 
4790= 
4800= 
4810=  £ 
4820=  1 
4830= 

484  0= 

485  0= 

486  0= 
4870= 
488  0= 
4390=  4 
4900=  3 
4910= 
4920= 
493 0=C 
494  0=C 
495 0=C 
496 0=C 
497 0=C 
498  0= 
4990= 
5000= 
5010= 
5020= 

5  03  0= 
504  0= 
5050=  £ 


PS I  FUNCTION 

SUFPOUT I HE  FPS I <M. NNMAX> 

COMPLEX  IX.F4.PSI 

COMMON  E 1 8  *  5 0  >  .  D  <  5  0 »  5 0  >  »  DNEG  <5  CO  .  P 1  C5  0>  .  F4  <5  0> 

1 .  SO  (  5  0)  .  F*  1 50!*  ?  SETA  * . 5  0 >  «  SET  AO  <50>  .  PS  I  *  5  CO1 

IX= < 0. .  1 . j 
MM=M+ 1 

PS  I  <MM>  =  < 0 .  .  o.:> 

DO  1  NN= 1 ?  NNMAX 
N=M+£^NN-£ 

PS  I  CMM>  =PSI  <MM>  + 1  OO  ♦SO  CNN)  ♦♦£^F4  <HN> 

CONTINUE 

RETURN 

END 

ASSOCIATED  LEGENDRE  POLYOMIALS.  OP.  OF  ARGUMENT  ETAO? 
orderpm.n;  with  n-m  even;  OF  TO  ORDER  N=M+£*NNMAX-£. 

SUBROUT I NE  POL YN  <ETAD , M . NNMRX> 

DIMENSION  PP  C3> 

COMPLEX  F4 

COMMON  E I G  <5  0>  >  D  <5 0 »  5 CO  .  DNEG  <5 CO  .  R 1  <5  CO  >  F4  <50> 

1 .  SO  <5 0 )  .  P  '-5 O'1 
SG'=SC!RT  <  1  -ETAD^ETAD) 

PP<1>=0 
PP  C£>  =  1 

IF  CM. EQ. 0>  GO  TO  1 
DO  2  L=1 . M 

PP  C£>  =  C2*L- 1 > ♦SQ^PP  C2> 

PCI)  =PP  <2> 

DO  3  NN=£ »  NNMRX 
N=M+2*NN-3 
DO  4  L=l.£ 

PP  C3>  =  C  C2+N- 1  >  ♦ETAO+PP  C2)  -  CN+M- 1  >  ♦PP  C 1  >  >  '  CN-M) 

N=N+ 1 

PP  < 1 >  =PP  <£> 
pp  <e>  =pp  <3> 

P  CNN) =PP C3> 

RETURN 

END 

DELATE  SPHEROIDAL  ANGULAR  FUNCTIONS.  S>  OF  ARGUMENTS 
C  AND  ETAO?  ORDER  M»N«  WITH  N-M  EVEN?  UP  TO  ORDER 
N = M + 2  ♦  N  N  M  A  X  -  2 . 

SUBPOUT I NE  OBANG CNNMAX . I RPMAX> 

COMPLEX  F4 

COMMON  E I G  <5  0>  ?  D  <5 0 »  5  0>  .  DNEG  <5  CO  .  R 1  <5  0)  t  F4  <5  CO 
1 9  SO  <5  CO  >  P  *5  ro  ,  SETA  <5  CO  » SET  AO  <5  CO 
DO  1  NN= 1 . NNMAX 
SETA  CNN)  =0 
DO  £  IRR=1 . IPPMAX 

SETA  <NN> =SETA CNN) +D  CNN. I PRO  ♦P<IRR> 
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50*0=  1  continue 

5070=  RETURN 

5080=  END 

5fiCM'i=C 

51 00=C  W  FUNCTION 

511 0=C 

5 1  £  0=  SUERDUT I NE  FW  CM  f  NNMAX> 

5 1 3 0=  COMPLEX  I X  f  F4 . W f  PS I 

514  0=  COMMON  E I G  <  5 0>  *  D  <  5  0  «  5  0>  « DNEG  <5  0)  «  P 1  <5  0>  «  F4  <5 O  ' 

5150=  SO  ('5 0>  f  R  <5 0>  *  SETR  <5 0>  f  SETRD  <5  CO  »  PS  I  <5  0>  >  W  *'5  CO 

5160-  IX=C0.fl.> 

5170=  MM=M+1 

5130=  WCMM>  =  <0.  f  0.) 

5190=  DO  1  NN= 1 » NNMRX 

5200=  N=M+£^NN-£ 

5£  1  0=  1  W  CMM>  =W  <MM>  ♦ 1  ^♦♦H^SETfiO  CNN>  ♦SD  CNN>  *F4  <NN> 

5££0=  RETURN 

5£3 0=  END 

5£4  0=C 

5£5 0=C  Y  FUNCTION 

5£*0=C 

5£7  0=  SUERDUT I NE  F Y  CM < NNMfiX> 

5£S 0=  COMPLEX  F4 , YO f YETRO fPSIt W 

5£9 0=  COMMON  EIG  <5 0>  » D  <5 0 >  5 0>  ?  DNEG  <5 CO  •  R 1  <5 CO  >  F4  <5 CO 

53  0  0=  1 9  SO  <5  CO  j  P  <5  CO  ?  SETR  <5  0>  •  SETfiO  (5  CO  >  PS  I  <5  CO  •  W  <5  CO  >  YC  5 

5310=  If  YETRO  <5  0> 

53£  0=  MM=M+ 1 

5330=  YQ<MM>«(0.  f  0.> 

534  0=  YETRO  CMM>  =  C  0 .  f  0 .  > 

5350=  DO  1  NN= If NNMRX 

5360=  N=M+£«-NN-£ 

537  0=  YO  <MM>  =YO  '•  MM>  +  C- 1  ♦♦N^P  1  <NN>  ♦SO  <NN>  ♦SETR  <NN>  +F4  CNN  • 

538  0=  1  YETRO  CMM>  =YETRO  CMM>  +  C- 1  '>  ♦♦N^R  1  CNN.)  ♦SETRD  CNN>  ♦SETR  >  * 

5390=  RETURN 

5400=  END 

541 0=C 

54£ 0=C  S  AND  U  FUNCTIONS 

543 0=C 

544  o=  SUERDUT I NE  FXU (MM *  MMMRX) 

545  0=  COMPLEX  I X  f  F4  f PS I • W f  YO f  YETRO  f  U ? X f  PT 

5460=  COMMON  EIG  <5 CO  f  D  (50 f  5 CO  «  DNEG  (50>  >  R1  <5 CO  f  F4  C5CO 

547  o=  1  f  SO  <5  CO  f  P  (5  CO  « SETR  <5  CO  f  SETRD  <5  0>  f  PS  I  <5  CO  f  W  (5  0>  f  YO  ■  5 

54g  o=  1  f  YETRO  (5  CO  f  U  (5  CO  f  X  (5  CO 

549  0=  I X=  <  0 . f 1 . > 

5500=  IF (MM. EQ. 1>  GO  TO  1 

5510=  IF  (MM. EC. MMMRX)  GO  TC  £ 

55£ 0=  PT=PS I (MM- 1 > ♦PS I (MM+ 1 > 

5530=  PTT=CRES  (PT> 

554  0=  IF  (PTT. EQ. 0. >  GO  TO  3 

555  0=  U CMM>  =£♦ I X ♦♦ (MM-£> ♦ (W (MM- 1 > +M (MM+ 1  >  >  / (PS I (MM- 1 > +PS T 

556  0=  X  (MM>  =£♦  I X^  (MM-£>  ♦  (W  (MM- 1  >  -l.l  (MM+ 1  >  >  /  (PS  I  (MM- 1  >  +PS  I 

5570=  RETURN 

553  0=  1  U  ( 1  >.=- 1  X^i.i  (£>  /PS  I  <£> 

5590=  X(l>  =  (0.f0.> 

5600=  RETURN 


5610=  2  U  (MM  >  =1V++  <MM-2)  ♦£♦1.1  (MM- 1  >  .'PS  I  <MM- 1  > 

562  0=  X  (MM) =U  (MM) 

5630=  RETURN 

56.40=  3  MMMAX=MM-  1 
5650=  RETURN 

5660=  END 

567 0=C 

568 0=C  Z  FUNCTIONS 

569  0=C 

57 0 0=  SUBROUT I NE  FZ  CMMMAX »Z»Zft»ZB» ZC » ZD) 

57 1 0=  COMPLEX  I X « 1 0 » Z  «. Zfi «  ZB n  ZC  n  ZD » F4 «  PS I • W n  YO  >  YET PD » U  t X 

572  0=  COMMON  E 1 6  <5  0>  »  D  <5 0-  50)  -  DNEG  <'5  0)  »R1  <5 0)  * F4  <5  0) 

5730=  1  f  SO  <5 0)  >  P  <50)  t  SETA  (50)  > SETRO  (50)  * PSI  (50)  » W  (50)  »  YO  (50) 

574  o=  1 ,  YETRO  (5 0)  j  IJ  (5  0)  * X  (5  0)  « CMPH I  (5  0)  »  SMPH I  (5  0) 

5750=  IX=C0.»1.) 

5760=  Z=  (0.  >  0. ) 

5770=  ZR=Z 

5730=  ZB=Z 

5790=  ZC=Z 

■so  fifi=  ■7ri="7 

5810=  MAXMM=MMMAX- 1 

5820=  DO  3  MM=1 n  MRXMM 

5830=  M=MM— 1 

584  0=  R=2 

5850=  E=1 

5860=  C=1 

5870=  IF  (MM. EQ. 1)  R=1 

5880=  IF  (MM.  EG!.  2)  B=£ 

589  0=  I F  (MM . EQ . 2)  C= 0 

5900=  IQ=IX^(-M> 

5910=  Z=Z+R^YETRO  (MM) ♦CMPH I (MM) 

5920=  IF  (MM. EQ. 1)  GO  TO  2 

593  0=  ZR=ZR+ 1 &♦ (U  (MM+ 1 ) ♦CMPH I  (MM+ 1 ) -B+U  (MM- 1 ) ♦CMPH I (MM- 1 ) ) 

5940=  1 ♦YO  (MM) 

595  0=  ZB=ZB+ 1 0^ (U (MM+ 1 ) ♦SMPH I  (MM+ 1 ) +U  (MM- 1 ) ♦SMPH I (MM- 1 ) ) 

5960=  1 ♦YD (MM) 

5970=  ZC=ZC+IQ^  (X  CMM+1)  ♦SMPHI  (MM+D-X  (MM-1)  ♦SMPHI  (MM-1)  ) 

5980=  1 ♦YO (MM) 

5990=  ZD=ZD+IQ+ (X (MM+1) ♦CMPH I CMM+1) +C^X  (MM-1) ♦CMPH I  (MM-1) ) 

6000=  l+YDCMM) 

6010=  GO  TO  1 

6  02  0=  2  ZR=Zfl+U (2) ♦CMPH I (2) ♦YO  ( 1 ) 

6  03  0=  ZB=ZB+U (£) ♦SMPH I (2) ♦YO  ( 1 ) 

6  04  0=  ZC=ZC+X (2) ♦SMPH I (2) ♦YD  C 1 ) 

6  05  0=  ZD=ZD+X (2) ♦CMPH I (2) ♦YD  ( 1 ) 

6060=  1  CONTINUE 

6070=  3  CONTINUE 

6080=  RETURN 

6090=  END 

6100=  SUBPOUT I NE  SCATTER ( VRP . ESCRTM . ANGLE* S I GMA) 

6110=  COMPLEX  A 1  *  P2 • B 1 » E2  *  H » P 1  * P2» R» U1 * U£ » X > Y » ESCAT*  J 

6120=  REAL  P I  f  THETA » An  Bn  D«  El n  E2 n  F 1 n  F2>  G1 « G2«  Cl n  C2*  KR* 

6130=  SS 1  n  S£  n  G1  n  T 1  n  T2  >  ARG «  V  *  W .  IERn  MMBS-J  0  n  MMBSJ 1 .«  MMBSJ2 

6140=  D I MENS I ON  VAR  ( 1 ) 

6150=  PI=3. 141593 
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1 

i  • 


1 


I 

I 


6160= 

THETR=  CVRR  <2>  >  -  360.  >  ♦  <2^PI> 

6170= 

kr=vrr  <■  1  > 

6180= 

R=l. 414214 

6190= 

B=—  <PI^4> 

62  0  0= 

!)=£♦  <  <  <2^Kfl>  /P I  >  ♦♦  C  1 .  /2 .  >  > 

6210= 

E1=CDS  <  CPI/-4.  >  -  <:THETp/£.  >  > 

6££0= 

E2=CGS  <  <P I  /A .- >  +  < THETP'  £ .  >  > 

623  0= 

F1=RBS<E1> 

624  0= 

F2=RBS  CE2!> 

6250= 

61=n^Fl 

6260= 

G2=D+F2 

6270= 

CALL  FRESNEL <G1» Cl > SI) 

6280= 

CRLL  FRESNEL  <62»  C£»  S£> 

6290= 

H=CMPLX  <0. »  B> 

6300= 

H=CEXP CH> 

6310= 

P1*GMPLX<C1»S1> 

632  0= 

P2=CMPLX<C£*S£> 

6330= 

fil=R^H*Pl 

6340= 

R2=R^H^P2 

6350= 

fi R6=£*Kfl^ <S I N <THETR> > 

6360= 

V=  a  +  SIN  (THETA)  >  ✓  <£■♦£  I N  <THETR>  > 

6370= 

1.1=  a  -  SI N  -::THETR>  >  /  <£♦£ I N  CTHETR)  > 

638  0= 

MMBS J2=  <£.-■  Rpi3>  ♦MMBSJ 1  CRRG » I  ER:<  -  MMBSJ  0  CRRG » I  ER> 

6390= 

X=  (  <R  1  ♦♦£•>  ♦W  +  CR2^2>  ♦V)  ♦MMBSJ  1  CRRG » I  ER> 

6400= 

Y=  <  CR  1  ♦♦£>  ♦l.l  -  <R2«"*2>  ♦V>  ♦MMBS J£ 

6410= 

J=  <  0 .  n  1 .  > 

6420= 

Y='v>J 

£.43  0= 

Y1=RERL CY> 

644  0= 

Y2=AIMAG<Y> 

6450= 

X1=RERL CX> 

6460= 

X2=RIMRG  <X> 

6470= 

ESCRT 1=X1-Y1 

6480= 

ESCRT£=X2-Y£ 

6490= 

ESCflT =CMPLX CESCRT 1 > ESCRT£> 

6500= 

ESCRT=ESCRT^J 

6510= 

ESCRT— ESCRT 

6520= 

ESCfiTF=REflL  CESCRT) 

6530= 

ESCRT I =R I MRG < ESCRT) 

654  0= 

ESCftTM  =  ccescrtr^2.)  +  cescrt  !♦♦£.))♦♦  a.  ✓£ 

.) 

6550= 

SIGMR  =  CESCRTM  ♦♦£.) 

6560= 

RNGLE=RTRN (ESCRT I ✓ESCRTR) 

6570= 

RNGLE  =  C  CRNGLE)  s  C2^F'I)  )  ♦BGO. 

6580= 

I F C CESCRTR . LT. 0. > . RND . CESCRT 1 . 6T . 0 . > )  RNGLE 

=RNGI 

6590= 

IF  C  <ESCRTR. LT. 0. )  . AND. CESCRT I. LT. 0.))  RNGLE 

=RNGI 

6600= 

RETURN 

6610= 

ENIi 

6620= 

SUBROUTINE  FRESNEL  CfiR,C*S) 

6630= 

Ii  I  MENS  I  ON  X  <2  0  0  0>  j  Y  <2  0  0 0>  •  ZX  <2  0  0  0)  •  ZY  C2 0 0  0> 

6640= 

INTEGER  NBIM 

6650= 

IFCCRR.GT.-O.  01)  .RND.  CRR.LT.  0.  01))60  TO  6 

6660= 

GO  TO  11 

6670=  6 

C=Q.  0 

6680= 

S=0.  0 

6690= 

GO  TO  4 

6700=  11 

H=.  01 

i 

I 


+  ISO.  I 
-  ISO.  ! 


I 

t 

] 

i 


I 
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6710= 

FP=hF3 • RR' 

!  6720= 

MI'!M=FI-'H  +  .1 

6720= 

K  =NDIN+1 

i  6740= 

DO  10  I  =  1  •  K 

j  675  0= 

H=I»h  -  H 

i  6760= 

X  I  ;>  =C03  1 . 57  08  *'i 

t  |  t*  i''  i*'  0= 

Y  <  I  >  =8 1 N  1 . 57  08  0^  <1 

6780= 

10 

continue 

6790= 

SUM£=0. 

i  6800= 

SIJM4=0. 

!  6810= 

1 

HH=.5»H 

!  682  0=C 

683  0=C 

INTEGRATION  LOOP 

'j  6840= 

HD  2  I=1*NPIM 

!  6850= 

SUM 1 =3UM2 

6860= 

SUM3=SUM4 

1 1  687  0= 

SUM2=SUM2  +  HH*<X< 

• '  688  0= 

SUM4=SUM4  +  HH^  <Y < 

h  6890= 

zxa>=suMi 

[i  6900= 

o 

c 

ZY  < I >  =SUM3 

f  j  6910= 

*t' 

ZX (NDIM> =SUM2 

K  6920= 

ZY <ND I M> =SUM4 

6930= 

7 

I F  <AFl .  GE  .0.1 )  GO  TO 

!  6940= 

C=  <- 1 . > ♦ZX  <NDIM> 

•  i  695 0= 

S=<-1.>»ZY<HDIM> 

|  6960= 

GO  TO  9 

[  6970= 

8 

C=ZX  <NDIM> 

6980= 

S=ZY <ND I M> 

6990= 

9 

CONTINUE 

7000= 

4 

RETURN 

7010= 

END 
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SCATTERING 

•HDBGE 


BY  Ft  METALLIC  CIRCULAR  BISK 
—  VERSION  1  £.■■" 1 7/73> 


<TYPE  KA=0  TO  STOP  PPDGPAM> 

<TYPE  KA=-1  FOP  A  BESPIPTICN  OF  THE  PARAMETERS* 
(NORMALIZATION:  ESCAT=A^EINC^ENORM/ <£^R> ♦EXP  <-J^K>R) > 

(ALL  ANGLES  IN  DEGREES* 


CM‘ STOPftGE 
2.1*8.  CP  SECONDS 
THETA  INCIDENT  =45 
POLRRIZRTION  =0 
THETR  SCRTTEPED  =45 
PHI  SCRTTERED  =0 


USED 

COMPILATION  TIME1 


WHICH  VARIABLE  IS  TO  BE  I NCREMENTEP? 1 
TYPE  NUMBER  OF  CASES:  £5 
WHAT  IS  THE  INCREMENT?! 


KA 

CROSS 

SECTION 

E 

NDPM 

SIGMA/  <PI^R^£* 

THETA 

PHI 

THETA 

PHI 

MAG 

PHASE 

MAG 

PHASE 

1. 00 

. 351E+00 

0. 

. 593E+00 

-££.65 

0. 

90.  00 

£.  00 

. 903E+0G 

0. 

. 950E+00 

-110.86 

0. 

90.  00 

3.  00 

. 263E+00 

0. 

. 513E+00 

148.85 

0. 

90.  00 

4.  00 

. 353E+00 

0. 

. 594E+00 

85.  18 

0. 

90.  00 

5.  00 

. 334E+00 

0. 

. 578E+00 

-3.41 

0. 

90.  00 

6.  00 

. 447E+00 

0. 

. 669E+00 

-105.68 

0. 

90.  00 

7.  00 

. 149E+00 

0. 

. 387E+00 

148.80 

0. 

90.  00 

8.  00 

. 738E-01 

0. 

. 272E+00 

58.53 

0. 

90.  00 

9.  00 

. 659E-Q2 

0. 

.81 £E- 0 1 

43.29 

0. 

90.  00 

10.  00 

. 194E-01 

0. 

. 139E+00 

-45.88 

0. 

90.  00 

11. 00 

. 534E-01 

0. 

. 231E+00 

-95.93 

0. 

90.  00 

1£.  00 

. 456E-01 

0. 

.£13E+00 

167 . 71 

0. 

90.  00 

13.  00 

. 10EE+00 

0. 

. 319E+00 

87.48 

0. 

90.  00 

14.  00 

. 627E-01 

0. 

. 250E+00 

-11.15 

0. 

90.  00 

15.  00 

. 982E-01 

0. 

. 313E+00 

-103. 23 

0. 

90.  00 

16.  00 

. 563E-01 

0. 

.237E+00 

1.94 

0. 

90.  00 

17.  00 

. 676E-01 

0. 

.260E+00 

94.77 

0. 

90.  00 

18.  00 

. 253E-01 

0. 

. 159E+00 

168.41 

0. 

90.  00 

19.  00 

. 275E-01 

0. 

.  166E+00 

-99. 14 

0. 

90.  00 

20.  00 

. 169E-01 

0. 

.  130E+00 

-51.91 

0. 

90.  00 

21. 00 

. 138E-01 

0. 

.  U7E+00 

45.27 

0. 

90.00 

22.  00 

. 294E-01 

0. 

. 171E+00 

97.  £2 

0. 

90.  00 

23.00 

. 199E-01 

0. 

. 141E+00 

-174.10 

0. 

90.  00 

24.  00 

. 513E-01 

0. 

. 226E+00 

-91.91 

0. 

90.  00 

25.00 

. 330E-01 

0. 

.  132E+00 

-5.  07 

0. 

90.  00 

Appendix  C 


Program  Listing  for  the  Fast  Fourier  Transform  (FFT) 

The  following  is  a  program  listing  of  the  Fast  Fourier 
Transform  which  performs  the  Direct  and  Inverse  Discrete  Fourier 
Transforms  (Ref  13:331). 

In  the  FFT  subroutine  argument  list,  X  represents  the  complex 
array  to  be  transformed,  M  represents  the  power  of  the  two  of  which  is 
the  length  of  the  FFT,  and  XI  indicates  whether  a  direct  FFT  or 
Inverse  FFT  is  to  be  performed.  If  XI  is  a  plus  one,  the  direct  FFT 
is  performed.  If  XI  is  minus  one,  the  inverse  FFT  is  performed. 

An  example  of  calling  the  FFT  subroutine  from  a  main  routine 
is: 

CALL  FFT  (DOG, 10,1.) 

In  this  Call  statement,  "DOG"  is  the  complex  array  to  be  transformed, 
"10"  is  the  power  of  two  of  which  is  the  length  of  the  FFT  and  "1." 
indicates  that  a  Direct  FFT  will  be  performed. 

A  few  precautions  In  using  the  FFT  subroutine  shown  should  be 
stated.  Using  the  above  example,  the  precautions  are: 

(1)  DOG  must  be  a  complex  array  with  2^  complex  locations. 

(2)  The  length  of  the  FFT  of  DOG  will  have  2^  complex  loca¬ 
tions. 


(3)  The  second  parameter  In  the  argument  list,  namely  “10" 
must  be  a  positive  integer  and  must  not  have  a  decimal 
point. 

(4)  The  third  and  last  parameter  in  the  list,  namely  "1."  must 
have  a  decimal  point. 

(5)  The  length  of  the  FFT  must  be  confined  to  half  the  maximum 
number  of  real  array  locations  that  are  possible  with  the 
computing  machine  used. 
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1 00= 

SUBPOUT I HE  FFT ‘X . M . X I > 

110= 

COMPLEX  XU>  •  U « to *  T 

120= 

130= 

HV£=H'-2 

140= 

NM1=H— 1 

150= 

J=1 

160= 

DO  7  1=1 » NM1 

170= 

IFCI.6E.J}  60  TO  5 

130= 

T=X < J> 

190= 

x<j:>=x<i> 

200= 

xa>=T 

210= 

5 

K=HV£ 

£20= 

6 

IFCK.GE.J)  60  TO  7 

£30= 

J=J-K 

£4  0= 

K=K/£ 

250= 

60  TO  6 

260= 

7 

J=J+K 

£70= 

PI =3. 14159265353979 

£30= 

DO  20  L= 1 ?  M 

£90= 

LE=2^L 

300= 

LE1=LE^2 

310= 

u=a.o»o.o> 

320= 

M=CEXP  CCMPLX  <  0 . » -X I  ♦P I  ''LE 1  >  > 

330= 

DO  £0  J=1*LE1 

34  0= 

DO  10  I=J*N-.LE 

350= 

IP=I+LE1 

360= 

T=X  < I P> ♦U 

370= 

X  <IP>  =X<D  -T 

330= 

10 

X  < I >  =X  <  D  +T 

390= 

20 

u=u*w 

400= 

IFCXI.6T. 0.>  RETURN 

410= 

DO  30  1  =  1,  N 

420= 

30 

X<I>=X<D/'N 

430= 

RETURN 

440= 

END 

Appendix  D 


Brief  Review  of  the  Selection  of  the  FFT  Parameters 

The  realization  of  the  Discrete  Fourier  Transform  on  the 
modern  digital  computer  in  the  form  of  the  Fast  Fourier  Transform  is  a 
valuable  tool  in  the  field  of  linear  systems.  However,  if  the  FFT 
parameters  are  not  carefully  picked  according  to  certain  rules,  the 
results  of  performing  an  FFT  on  a  sequence  will  provide  erroneous 
results.  So,  it  is  wise  to  explore  the  relationships  between  FFT 
parameters  in  the  time  and  frequency  domains  (Ref  15:282-284). 

First  of  all,  the  parameters  of  interest  are: 

T  -  increment  between  time  samples  (seconds), 
f  -  sampling  rate  (hertz)  s  1/T 
F  -  increment  between  frequency  samples  (hertz)  =  frequency 
resolution 

tp  -  record  length  (seconds)  =  effective  period  of  time  sig¬ 
nal  =  1/F 

fg  -  folding  frequency  =  f$/2  (hertz) 

fj,  *  highest  possible  frequency  in  spectrum  (hertz) 

N  -  number  of  samples  in  record 
In  order  to  avoid  aliasing. 


f.  >  2f. 


(D-ll 


For  a  desired  frequency  resolution 


tp  -  }  (0-3) 

If  f^  and  F  are  both  specified,  N  must  satisfy 


The  following  is  an  example  of  using  the  FFT  to  obtain  the 
DFT  of  the  106Hz  pulse  used  in  Chapter  III  and  Chapter  IV. 

The  pulse  consists  of  five  cycles  of  106Hz  R-f  centered  at 
the  origin.  The  pulse  is  shown  in  Figure  11.  Because  of  the  periodic 
assumption  of  the  time  signal  that  the  DFT  makes,  the  negative  half 
of  the  pulse  is  seen  to  be  shifted  up  in  time  ending  at  tp. 

The  first  step  is  to  estimate  what  the  highest  frequency  com¬ 
ponent  of  the  incident  pulse  is.  From  experiment  it  is  obvious  that 
the  Fourier  Transform  of  the  pulse  of  R-F  is  a  shifted  sinc(x) 
frequency  spectrum  which  never  dies  out  as  frequency  increases. 
However,  some  criteria  can  be  set  where  by  if  one  of  the  minor  lobes 
of  the  sinc(x)  waveform  goes  below  a  certain  percentage  of  the  main 
lobe  peak  value,  the  frequency  associated  with  that  minor  lobe  would 
be  considered  the  highest  possible  frequency  in  the  spectrum  of  the 


1 

pulse.  The  continuous  Fourier  Transform  of  the  10GHz  R-f  pulse 
is 


S 


(f) 


i-  v  in10  sin  h5.7  x  10~10(f-1010)1 
[15.7  x  10-10(f-1010)] 


(D-5) 


At  13GHz  up  from  the  center  frequency  of  10GHz,  or  23GHz,  the 
magnitude  of  the  sinc(x)  function  drops  to  5  percent  of  its  maxi¬ 
mum  value.  This  23GHz  will  thus  be  taken  as  the  highest  possible 
frequency  in  the  spectrum.  The  following  relationships  hold. 


fh  =  23  x  109Hz 
fs  >  2f h  =  46x1 09Hz 

T  <  0.0217  nsec 


(0-6) 


(0-7) 


If  f$  is  set  at  100GHz 

T  =  ■—  =  0.01  nsec  (D-9) 

Ts 

There  will  be  ten  samples  per  cycle  of  the  incident  pulse  and 
thus  50  samples  per  full  incident  pulse.  As  was  noted  in 
Appendix  C,  to  use  the  FFT  subroutine  shown  there,  the  length 
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of  the  FFT  must  be  a  power  of  2.  If  N  is  picked  at  1024: 


t  =  1024  x  T 

=  10.24  nsec 

(0-10) 

F  -  -r-  -  97.656  x  106Hz 

(D-ll) 

Therefore  the  highest  possible  frequency 

in  the  spectrum  of  the  pulse 

will  occur  at 

fh 

Sample  #  =  s  236 

(D-12) 

Therefore,  the  frequency  response  of  the  scattering  from  the  object 
the  incident  pulse  is  impinging  on,  can  be  cut  off  at  sample  250 
to  sample  1024-250=774  without  great  loss  in  accuracy.  Thus,  when 
the  DFT  of  the  pulse  and  the  sampled  frequency  response  are  multiplied 
and  inverse  transformed,  the  reflected  pulse  will  be  obtained. 

Because  of  the  way  the  points  were  picked,  the  process  of  multiplica¬ 
tion  of  transforms  and  inverse  transforming  will  perform  linear  con¬ 
volution  and  not  cyclical  convolution.  For  more  information  on  cylical 
convolution  refer  to  the  literature  (Ref  15:p. 284-295). 
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